Competing orders in a magnetic field: 
spin and charge order in the cuprate superconductors 
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We describe two-dimensional quantum spin fluctuations in a superconducting Abrikosov flux 
lattice induced by a magnetic field applied to a doped Mott insulator. Complete numerical solutions 
of a self-consistent large N theory provide detailed information on the phase diagram and on the 
spatial structure of the dynamic spin spectrum. Our results apply to phases with and without long- 
range spin density wave order and to the magnetic quantum critical point separating these phases. 
We discuss the relationship of our results to a number of recent neutron scattering measurements on 
the cuprate superconductors in the presence of an applied field. We compute the pinning of static 
charge order by the vortex cores in the 'spin gap' phase where the spin order remains dynamically 
fluctuating, and argue that these results apply to recent scanning tunnelling microscopy (STM) 
measurements. We show that with a single typical set of values for the coupling constants, our 
model describes the field dependence of the elastic neutron scattering intensities, the absence of 
satellite Bragg peaks associated with the vortex lattice in existing neutron scattering observations, 
and the spatial extent of charge order in STM observations. We mention implications of our theory 
for NMR experiments. We also present a theoretical discussion of more exotic states that can be 
built out of the spin and charge order parameters, including spin nematics and phases with 'exciton 
fractionalization '. 
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I. INTRODUCTION 

The determination of the ground state of the cuprate 
superconductors as a function of the hole density has 
been one of the central problems in condensed matter 
physics in the last decade. At zero hole density, it is 
well established that the ground state is a Mott insula- 
tor with long range magnetic Neel order. At moderate 
hole density, it is also widely accepted that the ground 
state is a G?-wave superconductor, all of whose important 
qualitative properties are identical those of the standard 
BCS-BdG theory. At issue are the ground states which 
interpolate between these well understood limits, and the 
manner in which they influence the anomalous properties 
at temperatures (T) above T c (the critical temperature 
for the onset of superconductivity). 

While a plethora of interesting proposals for these in- 
termediate states have been made, we will focus here on 
(in our view) the simplest possibility: the order param- 
eters characterizing the intermediate ground states are 
simply those of spin and charge density waves (SDW and 
CDW), and superconductivity (SC) itself. Apart from 
a small range at very low doping, which shall not be of 
interest in this paper, we know from neutron scattering 
experiments that there is SDW order collinearly polar- 
ized at the wavevectors 



(1.1) 



where a is square lattice spacing, and the wavevector shift 
from two sublattice order, < < 1/2, is a function of 
the doping concentration. In particular, strong motiva- 
tion for our study here was prmaded by the remarkable 
experiments of Wakimoto et alHu. They showed that the 
onset of superconductivity in L^-iSr^CuC^ occurs first 
at 8 = 0.055 (in a first-order insulator-to-superconductor 
transition) into a state which also has long-range spin- 
density-wave order at T = with a wavevector of the 
form ([□]) i.e. as 5 is moved away from the insulator 
at 5 = 0, the first conducting state is a SC+SDW state. 
As the ground state for large enough 8 is a SC state, it 
follows that there must be at least one quantum phase 
transition between the SC+SDW and SC states, and we 
will work with the simplest possibility that there is one 
direct transition at a critical 8 — S c . Wakimoto et al. also 
showed that such a transition associated with the vanish- 
ing of the SDW moment occurred for 8 C 0.14 (see Fig 
1 in Ref. g). We shall assume that the SC+SDW to 
SC quantum phase transition is second-order: direct evi- 
dence for critical magnetic fluctuations in L^-^Sr^CuC^ 
for 8 ~ 0.14 was provided in the neutron scattering ex- 
periments of Aeppli et aln. 

We will also discuss the appearance of local and long- 
range CDW order in the above phases. It is important 
to note that, throughout this paper, we use the term 
"charge density wave" (or "charge order") in its most 
general sense: such order implies that there is a periodic 
spatial modulation in all observables which are invari- 
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ant under spin rotations and time reversal, such as the 
electron kinetic energy, the exchange energy, or even the 
electron pairing amplitude; the modulation in the site 
charge density may well be unobservably small because 
of screening by the long-range Coulomb interactions. 

We note that the doping dependence of the magnetic 
order in the cuprates can be quite complex, varies signif- 
icantly between different compounds, and is influenced 
by the degree of disorder: the magnetic order may well 
be spin-glass like at the lowest energy scales at some 5. 
The SDW order is also enhanced in the vicinity of spe- 
cial commensurate values of the doping like 5=1/8 (see 
e.g. Fig 1 in Ref. ^), along with a suppression of SC 
order. In general, we do not wish to enter into most 
of these com plexities here, although we will mention (in 
Section I A) how our theory could be extended to explain 
the commensuration effects — som e other relevant issues 
will be discussed in Section VII. Our primary assump- 



tion is that the low energy collective excitations can be 
described using the theory of the vicinity of a quantum 
critical point between the SC+SDW and the SC phases; 
evidence supporting this assumption was also reviewed in 
Rcf. ^. This critical point is present either as a function 
of 5 in the material under consideration, or in a gener- 
alized parameter space but quite close to the physical 
axis. 

It is also important not to confuse this magnetic quan- 
tum critical point, with other proposals for quantum crit- 
ical points near ontimal doping that have appeared in the 
recent literatureM. These latter critical points are near 
5 « 0.19, and are probably not associated with long- 
range spin density wave order at a wavevector of the 
form (Q). This paper will discuss magnetic transitions 
at smaller doping. 

Upon accepting the existence of a second order quan- 
tum critical point at T = between the SC+SDW and 
SC phases, a powerful theoretical tool for the analysis 
of experiments becomes availablcQ. The structure of the 
critical theory, and its associated classification of eigen- 
perturbations, allows a systematic and controlled theory 
of the spin excitations in the SC and SDW phases on 
either side of the critical point. Such an approach was 
recently exploited to study the influence of non-magnetic 
Zn and Li impurities in the SC phasea. In this paper we 
will use the same tools to study the influence of an ap- 
plied magnetic field, oriented perpendicular to the CuC>2 
layers, on both the SC and the SC+SDW phases. An 
outline of our results, has already appeared in previous 
communicationsailj'EJI: here we will present the full nu- 
merical solution of the our self-consistent equations for 
the dynamic spin spectrum in an applied field, along with 
a number of new results. Measurements of the spin and 
charge correlations in the presence of such an applied 
magnetic field have appeared r ecently in ajiumber of 
illuminating neutceii scattering,E3t§EJ NMRMEj'EZI and 
STM experimentsta, and we will compare their results 
with our prior predictions. 



A. Order parameters and field theory 

The field theory for a SC to SC+SDW transition in 
zero applied magnetic field can be expressed entirely in 
terms of the SDW order parameter which we will intro- 
duce in this subsection; the quantum fluctuations of the 
SC order can be safel y ne glected, a point we will dis- 
cuss further in Section VII . Consideration of the applied 



magnetic field will appear in the following subsection. 

We introduced above the wavevectors of the SDW or- 
dering K sx and K S j,; almost all of our analysis will apply 
for general values of but the value d = 1/8 is of par- 
ticular interest above a doping of about 1/8. To obtain 
an order parameter for such a SDW, we write the spin 
operator S a (r, r), a = x, y, z, at the lattice site r as 



S a (r,r) = Re[e iK -- r $ xa (r, 



(1-2) 

where <& x ,ya are the required order parameters. Except 
for the case of two sublattice order with $ = (which 
we exclude for now), the fields <& x . ya are complex. These 
fields can describe a wide variety of SDW configurations, 
but we now list the two important limiting cases. 
(i) Collinearly polarized SDWs, for which 



<Mr,r) = e w < r ' T >n a (r J T) 



(1.3) 



where n a is a real vector and is also real (and similarly 
for & xa ). Parameterized in this manner, and for r? a = 
constant (summation over the repeated index a is implied 
here and henceforth), the order parameter Q ya belongs 
to the space (5*2 x S\)/Z2, where S n is the n-dimensional 
surface of a sphere in n + 1 dimensions, and Z v is the 
discrete cyclic group of p elements. The Z2 quotient is 
necessary because a shift 9 — ► 8 + it is equivalent to a 
rotation which sends n a — > —n a . 
(ii) Circular spiral SDWs, for which 

$j, a (r,r) = n la (r, r) +m 2Q (r,r) (1.4) 

where n\^a are t w0 rea ^ vectors obeying n\ a = n^ a 
and n\ a nia = (and similarly for <5> X a)- Now for 
n\ a = constant, the order parameter $> ya belongs to the 
space SO(3)^ S 3 /Z 2 (see ^.-Section 13.3.2 in Ref. 
The experimental evidencelial23 supports the conclusion 
the SDW ordering in the cuprates in collinear, but the 
present formalism allows a common treatment of both 
the collinear and spiral cases. This complex- vector for- 
mulation of the SDW order allows treatment of the SDW 
quantum transition by a straightforward generalization 
of the real-vector theory used for the Neel state in the 
insulator; r *elated pointsiave been made by Castro Neto 
and HoneEil and ZaanenE3. The same approach was also 
used by Zachar et alr3 to treat the onset of SDW order 
at finite temperatures, as we will indicate below. 

Along with the SDW order, CDW order may also ap- 
pear. We parameterize the charge density modulation 
by: 

5p(r, r) = Re [e iK -' r ^(r, r) + e iK ^ r y (r, r)] (1.5) 
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where ¥L cx . y are the CDW ordering wavevectors and 4> x ,y 
the corresponding complex order parameters. The quan- 
tum numbers of the obse rvab le Sp are identical to those of 
S^, and so by squaring ( |l p2^ we see that associated with 
the SDW is a CDW withB K cx = 2K SX , K cy = 2K 
(modulo reciprocal lattice vectors), 
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<fe(r, t) oc ^(r, r), and (f> y (r, r) oc $> 2 ya (r, r). (1.6) 

Note that this CDW is absent for the case of a circular 
spiral SDW (in which case <E>^ ya — 0) but is necessarily 
present for a collincar SDW. In principle, in a state with 
condensates of both $> xa and & ya , a CDW can also be 
present at wavevector K SX +K SJ/ ; we will not consider this 
possibility here as it does not seem to be experimentally 
relevant. As was emphasized in the third paragraph of 
Section |, we are using the term CDW here in its broadest 
sense: there is a modulation at the wavevector K c in all 
observables which are invariant under spin rotations and 
time reversal. The precise nature of the CDW order may 
be determined from an analysis of the STM spectrum- 
this has been discussed recently in Refs. |24| , p5| . 

The order parameters & x , ya , <fix,y allow a rich variety 
of phases and phase transitions in the presence of back- 
ground SC order. These will be discussed in some detail 
in Section VI. Central to a description of these phases 



is an understanding of the symmetries respected by any 
effective action for the order parameters. We describe 
these below and then focus on a particular phase transi- 
tion of physical interest. 

An obvious symmetry is that under spin rotations; this 
is described by the group SU(2), and the fields $ Xj2/a 
transform as <S* = 1 vectors labeled by the index a. In 
addition, there is an independent sliding symmetry 



3> 



x,ya 



x.ya 



(1.7) 



associated with the translational symmetry of the un- 
derlying lat tice model: t rans lating r to r + (ma, 0) (m 
integer) in (L2) leads to ( 1.7 ) wi th 9 X = rmr(l — 2$) and 
y = rmr (i? was defined in (1.1)). For $ irrational, we 
see that all real values of 6 x<y can be generated with the 
different choices for m, and hence the sliding symmetry 
is U(l)xU(l). For rational with 1/2 - i? = p'/p, and 
p' , p relatively prime intege rs, only integer multiples of 
9 x ^ y = 2ir/p are allowed in (1.7); in this case the sliding 



symmetry is reduced to Z p x Z p . The difference between 
U(l) and Z p will not be material to any of our results for 
p > 2. In a similar manner, we can also determine the 
action of other elements of the square lattice space group 
on <& x ,ya and we mention two important cases: under a 
spatial inversion we have <fr x ,ya —> & x ya i and under the 
interchange of x and y axes, we have $ m <-> $> ya . 

We now apply these symmetries to determine the ef- 
fective action of a physically relevant transition discussed 
earlier in the introduction (and in the phase diagrams of 
Section ||): that between the SC+SDW and SC phases. 
This transition is driven by the condensation of $> x . va ] 
if the SDW order is collinear, it will drive a concomi- 
tant CDW order, as discussed above. Supplementing the 



symmetries by a renormalization group (RG) procedure 
which selects terms with smaller powers of <&rx-n 
fewer spatial and temporal gradients, we obtainlilJr 1 ' ~ '' 
the effective action 



<S<i> 



\d T ^ xa \ + w? \d x $> xa \ + v\ |c\$ a 



d 2 rdr 



Vl \d x $y 



Vf\dy$ ya \ 



\* 2 ya \ )+w 1 \$ xa \ 2 W 



(1.8) 



Note that first-order temporal gradient terms 
like $>* xa d-r& xa are forbidden by spatial inversion 
symmetryES. In principle, first-order spatial gradient 
terms like i$* xa d x Q xa are permitted by all symmetries; 
such terms lead to a shift in the wavevector at which 
SDW fluctuations are largest, and we assume that they 
have already been absorbed by our choice of K sx . Here 
vi and V2 are velocities, which are expected to be of 
order the spin- wave velocity, v, of the Neel state in the 
undoped insulator. The parameter s tunes the system 
from the SC phase (s > s c ) to the SC+SDW phase 
(s < s c ), where s = s c is the non- universal location 
of the quantum critical point between these phases; 
experimentally, s can be varied by changing the doping 
concentration. The action also contains a number of 
quartic non-linearities: the RG analysis shows that these 
are strongly relevant perturbations about the Gaussian 
theory, and will play a crucial role in our analysis below. 
The coupling ui selects between the collinear and spiral 
SDW states: for 112 > 0, the circular spiral state (which 
has = 0) is selected, while 112 < prefers a collinear 
SDW. The couplings wi.2,3 lead to correlations between 
the orders at K sx and K SJ( — if these are attractive, 
the s < s c phase will have simultaneous orderings 
at both wavevectors, and spatial pattern will have a 
checkerboard structure. 

We have also neglected the couplings to the low en- 
ergy nodal quasiparticles, which are additional excita- 
tions of the SC phase carrying spin; their effects are sup- 
pressed by the constraints of momentum conservation, 
as they can damp the $ quanta effectively only if K SIja 
equal the separation between any two nodal points. The 
case where this nesting condition is satisfied has been 
considered earlierc3, but we will not enter into it here 
for simplicity: essentially all of our results here on the 
phase diagram in an applied magnetic field apply also 
to the case where the nesting condition is obeyed. For 
completeness, in Appendix H we also discuss the role of 
spin symmetry breaking Dzyaloshinskii-Moriya interac- 
tion present in La^-aSraCuO^. We show that it helps 
stabilize collinear SDW order in a certain direction; how- 
ever its effect is very small and will be neglected in the 
rest of this paper. 

For the particular rational value $ = 1/8, the 
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U(l)xU(l) sliding symmetry is reduced to a discrete 
Z s x Z 8 symmetry under which 9 x y in (1.7) are only 

allowed to be multiples of 7r/4. This reduced symme- 

whosc structure has 



(we use units with h = fee = 1 throughout) 



1. 



try allows additional term s in (] 

been discussed earliercllEii Such terms help-choose be- 
tween site- and bond-centered density wavesEJ, and could 
also lead to the enhancement of the moment observed by 
Wakimoto et aln near S — 1/8. However, these terms 
are very high order (eighth) in the $ fields, and conse- 
quently they have a negligible effect on the issues we are 
interested in here: so we will not consider them further. 

It is useful to compare our treatment here of the 
SC+SDW to SC transition with others in the literature. 
It is essential for our purposes that the spin/charge order- 
ing is taking place in a background of SC order, as that 
gaps out the fermionic excitations except . p o ss ibly at spe- 
cial points in the Brillouin zone. TheoriesE§E3 which con- 
sider SDW/CDW order in a Fermi liquid have additional 
damping terms in their effective action which change the 
universality class of the transition, change the dynamic 
exponent to z = 2, and do not obey strong hyperscal- 
ing properties as the quartic couplings are marginally 
irrelevant in this case. We have also taken a genuinely 
two-dimensional view on the SDW/CDW ( "stripe" ^-fluc- 
tuations in our approach. An alternative approaches as- 
sumes there are intermediate scales on which the physics 
of the one-dimensional electron gas applies, although a 
crossover to similaii_two-dimensional physics occurs on 
large enough scales^. 



B. Influence of an applied magnetic field 



An applied magnetic field has a Zeeman coupling to the 
spin of the electrons, and this is present for any direction 
of the applied field. However, the Zeeman splitting of 
the magnetic levels has only a minor effect, and can be 
safely neglected compared to the much stronger effects 
near s = s c that we consider below. We discuss the 
influence of the Zeeman term in Appendix |b], and will 
not consider it further in this paper. 

The dominant effect of the field is via its coupling to 
the orbital motion of the electrons, which is sensitive only 
to the component of the field orthogonal to the layers. 
The reason for this strong effect is simple: there is SC 
order in the orbital wavefunction of the electrons, and 
the diamagnetic susceptibility of the SC state to the ap- 
plied field is infinite. However, as the SC order is non- 
critical across the transition at s = s c , it is mainly a 
quiescent spectator and its response can justifiably be 
treated in a static, mean-field theory. Consequently, we 
model the complex SC order parameter tp(r) in the fa- 
miliar Abrikosov theory with the free energy per layer 



T= I d 2 r 



- a|V(r)| 2 + ^(r)| 4 



2m* 



A V(r) 



.(1.9) 



Note that unlike <5> x ,ya, tp is not a fluctuating variable, 
and described completely by its mean-value (which will 
be r dependent). We will work entirely in the limit of ex- 
treme type-II superconductivity (with Ginzburg-Landau 
parameter kql = oo); so there is no screening of the mag- 
netic field by the Meissner currents, and V r x A = Hz, 
the applied, space-independent magnetic field. 

To complete the description of the model studied in 
this paper, we now need to couple the SC and SDW or- 
der parameters together. The simplest term allowed by 
symmetry is a connection between the local modulus of 
the order parameters: 



S<S>tp — K 



d 2 rdr|V(r)| 2 (|^ a (r, t)\ 2 



I* 



ya(r,T)\' 



(1.10) 

For k > 0, we can induce a competition between the 
SC and SDW orders, in that the SDW order will be en- 
hanced where the SC order is suppressed and vice-versa. 
The microscopic origin of the coupling k is discussed in 
Appendix ^|. 

Although 5$^ will be the primary coupling between 
the SDW and SC orders, an additional allowed term will 
be important for some purposesEU. To understand this, 
notice that all terms in 5$ and 5$^, are invariant under 
the sliding symmetry ([l]?]). This means that, with the 
present terms, the CDW order is free to slide arbitrarily 
with respect to any vortex lattice that may be present 
in the SC order tp. This clearly cannot be true, as lat- 
tice scale effects should pin the two modulations with 
respect to each other. The simplest additional coupling 
which will provide this pinning can be deduced by notic- 
ing that there shou ld be a coupling between the charge 
modulation in (1.5) and the local modulus o£-the super- 
conducting order; this is induced by the termEil 

Slat = -C£ / rfr^(r)| 2 Re[e lK - r $L(r,r) 

+ e *-%(r,r)]. (1.11) 

Notice that we are now performing a discrete summation 
over the lattice sites r, rather than integrating over a spa- 
tial continuum: this is a direct consequence of the rapidly 
oscillating factors e lKcx ' r and e jKc »' r which do not have 
a smooth continuum limit. Indeed, in regions where ip( r ) 
is smoothly varying, these rapidly oscillating factor will 
cause the summation over r to vanish. So the expres- 



sion (1.11) is appreciable only over regions where 
is rapidly varying, and this happens only in the cores of 
the vortices. As the centers of the vortices are identi- 
fied by the zeros of ip(r), and we are mainly interested in 



scales larger than vortex core size, we can replace (1.11) 
by the following expression, which is-.morc amenable to 
an analysis in the continuum theorytJ: 



Si 



lat 



C J drRe[e^ (t>l a (r v ,T) 



r„,i/>(r„)=0 



+<f> 2 ya {r v ,r))] 



(1.12) 



Here the summation is over all points r„ at which ip{v v ) = 
(these are the centers of the vortices), and w is a phase 
which depends upon the microscopic structure of the vor- 
tex core on the lattice scale. The action S\ a t is not in- 
variant under the sliding symmetry, and so will pin the 
CDW order. 

We are now in a position to succinctly state the field- 
theoretic problem which will be addressed in this pa- 
per. We are interested in the partition function for 
SDW/CDW fluctuations defined by 



Z[1>(t)] = / D$ M (r,T)0$ 9a (r,r) 



x <'-M> ( --= - <S* - - <Sj at )(1.13) 



accompanied by the solution of 
(51nZ[^(r)] 



(1.14) 



which minimizes — InZ [V>( r )] to determine the optimum 
ip(r). Note the highly asymmetrical treatment of the 
SDW and SC orders: we include full quantum- mechanical 
fluctuations of the former, while the latter is static and 
non-fluctuating. This asymmetry is essentially imposed 
on us by the perspective of magnetic quantum critical- 
ity, and the fact that we are developing a theory of the 
SC+SDW to SC transition. This asymmetry should also 
be contrasted with the symmetric treatment o£-SC and 
SDW quantum fluctuations in other approaches^ 2 ! 



C. Physical discussion 

The primary purpose of this paper is to determine the 
phase diagram and low-energy spectrum of SDW and 
CDW fluctuations of Z as a function of the applied .field. 
H. A summary of our results has already appearecfJlillO 
and detailed numerical solutions appear in the body of 
the paper; here we expand on the central physical idea 
to provide an intuitive understanding of our results to 
readers who do not wish to study the details in the re- 
mainder of the paper. We will initially ignore the pin- 
ning described by <5>i a t, but will discuss its consequences 



in Section I C 1 



Let us begin in the SC phase with s > s c and con- 
sider the $ XQ fluctuations in a simple Gaussian theory 
(the considerations of this subsection apply equally to 
$ yQ , which we will not mention further). Assume ip{ T ) 



(Energy) 
A 




FIG. 1: A sketch of the potential Vo(r) (thick full line) in 
the presence of a vortex lattice. Also shown is the exciton 
wavefunction Ho(r) which solves (1.17) for V(r) = Vo(r) with 
eigenvalue A 2 . Note that there is no drastic change in this pic- 
ture as A 2 \ 0: the peaks in Ho(r) remain exponentially lo- 
calized within each vortex core, on a length scale much smaller 
than the vortex lattice spacing. We argue in the text that 
strong interaction corrections to Vo(r) invalidate this form 
for So(r) and the correct structure is shown in Fig. 0. 



has been determined by the minimization of J 7 , and so 
takes the standard form in an Abrikosov flux lattice. The 
Gaussian fluctuations of $ M are described by the effec- 
tive action 



Sr 



d 2 rdr 



Vl\d x $ xa f +vl\dy$ xa \ X 



+ V(r)|* I 



(1.15) 



To leading order, the effective potential V(r) is given by 
V = Vo where 



V (r) = s + «#(r)| 2 



(1.16) 



A sketch of the spatial structure of Vb(r) is shown in 
Fig [|: because ip(r) vanishes at the centers of the vor- 
tices, Vrj(r) has well-developed minima at each such 
point. Indeed, there can even be regionSrin each vor- 
tex care where Vo(r) < 0, and Arovas et alS2 and Bruus 
et alEj argued that superconductivity would 'rotate' or 
transform into static Neel oijder in such a region. In our 
treatment of dynamic SDWcH, we see that the struc- 
ture of the magnetism is-determined by the solution of 
the Schrodinger equatior£j 



{-vldt 



^ 2 + V(r))~ (r) = A 2 ~ (r) (1.17) 



where Srj(r) is the lowest eigenmode of ( 1.17 ), the eigen- 
value A 2 is required to be positive for the stability of 
the Gaussian theory Sg- The energy A is the spin-gap, 
and So (r) then specifies the envelope of the lowest energy 
SDW fluctuations; in other words Srj(r) is the wavefunc- 
tion of a S = 1 exciton associated with dynamic SDW 
fluctuations. Note that A 2 can be positive even if there 
are regions where V(r) < 0. A sketch of the spatial form 
of So(r) is shown in Fig [l] for a particular small value of 
A 2 and V(r) = Vrj(r). Observe that Sn( r ) is peaked at 
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the vortex centers, but decays rapidly outside the vortex 
cores over a SDW localization length £ ~ Vi^/VVl — A 2 , 
where Vt is the value of Vo(r) outside the vortex cores 
(see Fig H). 

Remaining within the Gaussian theory specified by 
(1.15) and (1.16), we now consider the consequences 



(Energy) 

A 



of raising the value of H in the hope of reaching the 
SC+SDW phase. With increasing H, the vortex cores 
will approach each other, and we expect that the value of 
A 2 will decrease. Indeed, the picture of Fig |l| holds all the 
way up to the point A = 0; beyond this field the Gaus- 
sian theory becomes unstable and this signals the onset of 
the SC+SDW phase driven by the condensation of Q X a- 
Note that the localization length £ ~ V12/VVL — A 2 of 
the SDW order peaked in the vortex cores remains finite 
all the way up to the critical point. This localization 
length £ must be clearly distinguished from the spin cor- 
relation length, £ s : the latter is associated with correla- 
tions between different vortices, and arises because there 
is an exponentially small coupling between magnetism 
in neighboring cores. Thus this simple Gaussian theory 
yields a picture of dynamic magnetism appearing first in 
the vortex cores, with possible weak correlations between 
neighboring cores. Such a viewpoint was also discussed 
by Lake et al£3 who proposed "spins in the vortices" but 
noted that the large value of £ s implied coupling between, 
nearby vortices. Following our worko, Hu and ZhangLj 
also presented a picture of dynamic SDW fluctuations 
similar to the one above. 

We now argue that corrections beyond the Gaus- 
sian theory approximation invalidate the above picture 
when A becomes smalltl. Indeed, the picture of nearly- 
independent, localized magnetic excitations in each vor- 
tex core holds only then A is of order the spin exchange 
energy J; such high energy magnetic excitations are ex- 
pected to strongly damped by the fermionic quasiparti- 
cles. Also, the validity of the present continuum model 
is questionable at scales as short as vortex core size and 
at energies of order J: a full solution of the BCS theory 
of the underlying electrons is surely needed, and sub- 
sidiary order parameters may well develop within the 
vortex cores. However, as A is lowered, we will now 
argue that the physics is actually dominated by the large 
region outside the vortex cores, where the present contin- 
uum approach can be used without fear, and the subtle 
issues of the short-distance physics within the core can 
be sidestepped. The central weakness in the analysis of 
the previous paragraph is that it does not account for the 
repulsive interactions between the bosonic <fr xa ex- 
citon modes that are condensing .As has been discussed 
in different contexts long agotSEa, such interactions are 
crucial in determining the structure of the lowest energy 
state in which condensation occurs. In particular, it is 
well known that the effect of interactions is to delocalize 
the lowest energy states: bosons initially prefer to oc- 
cupy strongly localized, low energy states, but then their 
repulsive interaction with subsequent bosons drives the 
energy of such states up. Bray and MooreEa presented an 




FIG. 2: A sketch of the potential Vo(r) (thick full line) in 
the presence of a vortex lattice along with t he tru e form of 
the exciton wavefunct ion 5 p(r) which solves (1.17) with the 
full potential V(r) in (1.18). The spatial structure of Ho(r) as 



A \ is characterized by the vortex lattice spacing. 



argument demonstrating that in the vicinity of the con- 
densation, the localization length must diverge as one 
approached the bottom of the band of states of interact- 
ing bosons in the presence of an external potential. To 
apply th eir a rgument in the present context, we need to 
replace (1.16) by 



So 



V(r) = V (r) + (l^±^)(|^(r,.)| 2 ) 

= S + *#(r)| 2 + (4U1 + 2U2) (|$, Q (r,r)| 2 ) 5G ; (1.18) 

the additional terms arise from a Hartree-Fock decou- 
pling of the quartic interaction terms in <S$, and the ex- 
pectation values have to be evaluated self-consistently 
under the Gaussian action in (1.15) which itself depends 
upon V(r). Note that the perspective of magnetic criti- 
cality requires that we account for the u% t 2 interactions, 
as these are strongly relevant perturbations about the 
Gaussian theory; so we are led to ( 1.18| ) also by a naive 
application of the RG approach. We will present de tailed 
numerical solutions of equations closely related to (ITS) 
in the body of the-pape r. An adaption of the argument of 
Bray and Mooret3 to ( 1.18 ) was given in Ref . ||, and we 
will not repeat it here: the main result is that the length 
scale £ characterizing the lowest energy state So(r) can- 
not remain finite as A \ 0. Instead the states around 
neighboring vortex cores overlap strongly, and Sn(r) is 
characterized by the vortex spacing itself. A sketch of 
the actual structure of S ( r ) is shown in Fig g. The spin 
correlation length, £ s , does not have a direct connection 
with the spatial form of Sn(r) itself, but is instea d rela ted 
to an integral over a band of states which solve ( 1.17 ) at 
finite momentum, as we shall discuss in Section I C 1 and 
later in the paper. 



It i s worth noting here that the passage from (Lit) to 
( 1.18 ) in zero field is precisely that needed to reproduce 
the known properties of magnetic quantu m cr itical points 
in other situations. In one dimension, (1.16) would im- 
ply th at the re is no barrier to magnetic long-range order, 
while ( 1.18| ) correctly implies that the presence of the 
Haldane gap, and reproduces its magnitude in the semi- 
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classical limitfl. At finite temperature, ( 1.18p yields the 
correct crossovers in the magnetic correlation length in 
the vicinity of the spin ordering transition in two dimen- 
sions. Although w e will not present detailed solutions 
on this case here, (1.18) is also expected to provide a 
reasonable description of the magnetic crossovers at fi- 
nite temperatures in the vicinity of the SC+SDW to SC 
transition in the presence of a magnetic field. 

With the knowledge of the spatial structure of the ex- 
citon wavefunction S (r) in Fig ||, the origin of our main 
resultaJ can be easily understood. As the vortex cores oc- 
cupy only a small fraction of the system volume, the mag- 
nitude of the energy A 2 is influenced mainly by the struc- 
ture of ip(r) in the remaining space. Here, the predom- 
inant consequence of the magnetic field is the presence 
of a superflow with velocity v s = —5T/8A. circulating 
around each vortex core. Focusing on the region around 
a single vortex at the origin r = (0,0), the superflow 
obeys |v s | ~ 1/r in the wide region £ < r < L„ where 
£o = l/ V2m*a is the vortex core size, L v ^ 



¥ H/c) 



-1/2. 



so the average superflow kinetic energy is 



r 


d 2 r 


J to 




f 


d 2 r 







where H® 2 is the upper critical field for the destruction of 
the Meissner state at the coupling constant correspond- 
ing to the point M in Fig [| below. This kinetic energy 
is a scalar with the same quantum numbers and symme- 
try properties as |?/>| 2 : henc e, via the coupling in Sg,^, 
in (1.10), the value of ( 1.19| ) feeds into all the effective 
coupling constants in 5$ in (|k^). The most important 
modification is that the tuning parameter s gets replaced 
by 



H / M 
Scff( #) = s _C_ln(^ 



(1.20) 



where C is a constant of order unity. The imp licati on of 

( 



1.20 ) is that we may as well replace V(r) in ( 1.15 ) and 



1.18) by 



V(r) « s eff (ff) 



(1.21) 



to obtain a first estimate of the consequence of the mag- 
netic field in the vicinity of the S C+S DW t o SC transi- 
tion. The H dependence in ( [1^0 ) and (1.21) is sufficient 
to determine our main resultstil: the small H portion of 
the phase diagram in Fig |[ the intensity of the elastic 
scattering Bragg peak in the SC+SDW phase, and the 
energy of the lowest energy SDW fluctuation in th e SC 
phase. In particular, it follows directly from ( 1.20] ) that 
the small H portion of the AM phase boundary in Fig || 
between the SC and SC+SDW phases behaves as 



H 



2(s-a c ) 
k1u(1/(s - s c )) ' 



(1.22) 



H 



"Normal" 
(Charge order) 
D 




FIG. 3: Zero temperature phase diagram as a function of the 
coupling s and the magnetic field H in an extreme type II 
superconductor described by (1.13). The theory is accurate 
in the region of small H, and only qualitati vely correct else- 
where (H is measured in units described in ( |3.2| ) . The phases 
without SC order are likely to be insulators, and the "Nor- 
mal" phase is expected to have residual CDW order, which 
is initially induc ed by the pinning terms in 5i at as discussed 
in Section J C % The positions of the phase boundaries are 

The path Pi denotes the loca- 



summarized in Section 



III 



tion etthe original neutron scattering measurements of Lake 
et al£3, and the path P2 the subsequent neutron scattering 
measurements of Khaykovich et alx3 amd Lake et al&3. The 
STM measurements of Hoffman et al.tB are also along path 
Pi. 



Note that this phase boundary approaches the s = s c , 
H = quantum critical point with vanishing slope. This 
implies that a relatively small H for s > s c will suc- 
cessfully move the system close the AM phase boundary, 
and so produce low energy spin excitations. This should 
be contrasted to the corresponding ii-dependent phase 
boundary of the SDW phase in insulators which is dis- 
cussed in Appendix here, there is no orbital diamag- 
netism, only the Zeeman coupling is operative, and the 
phase boundary approaches H = with infinite slope. 
Evidently, Zeeman effects are much weaker and can be 
justifiably neglected. 

We conclude rH^r^i^^iijmJj^^ of 
earlier workSfflS^ffiHSWffiH on vortex 
magnetism and the change in perspective that has been 
offered h 
Sachdevt 



are by our analysis. 



was proposed in by 



and Nagaosa and Leeci that vortex cores in 
the underdoped cuprates should have spin gapj-corre- 
lations characteristic of Mott insulators. ZhangLJ and 
Arovas et alEB described vortex core correlations in terms 
of static Neel order, and estimated that the field-induced 
moment would be proportional to H (in our phase di- 
agram in Fig ||, vortices with static moments are only 
present in the S C+S DW phase, and as we will revie'st 



below in Section VP , the average moment increases asEl 
H \u(\/H) for small H). Our discussion here also uses 
the SDW order parameter, but allows it to fluctuate dy- 
namically into a spin gap state, and so interpolates be- 
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tween these earlier works. A separate description of ,vor- 
tex cores in terms of 'staggered flux' correlationsCJO has 
also been proposed. One of our central points here is 
that while the vortex core correlations may well be quite 
complicated (they are dependent on lattice scale effects, 
and difficult to distinguish from each other as the short- 
distance 'order' fluctuates dynamically), these issues can 
be sidestepped: a reliable continuum theory can be de- 
veloped by considering first the dominant effects arising 
from the interplay between superconductivity and mag- 
netism in the superflow region outside the vortex cores. 
Spin density wave correlations induced in these regions 
may leak into the vortex cores, but our treatment is not 
expected to be reliable in the latter region: the nature of 
the electronic correlations in the vortex cores remains an 
open question. 

Our continuum treatment of dynamic and static spin 
density wave order differs from earlier works in several 
key aspects. An important feature of Refs. |3^|34| is the 
static mean-field treatment of the SDW order in the vor- 
tex cores, which is imposed by their "SO(5)" picture of 
SC order outside the core&,'rotating' into static antifcr- 
romagnetism in the coresE3. This should be contrasted 
to our approach, in which magnetic quantum criticality 
implies dynamic magnetic fluctuations while the SC or- 
der can be safely considered static. Further, Refs. |33j , |34| 
assumed the (near) equality of the gradient and "mass" 
terms for the SC and two-sublattice SDW order parame- 
ters, as naturally suggested by the dynamic SO (5) sym- 
metry, which requires a symmetry between the excited 
states in the SC and SDW phases. As a result they 
found static two-sublattice magnetization induced by the 
vortex core, over a scale which was of order the vor- 
tex core size, £o, an d in a regime where superconduc- 
tivity was essentially completely suppressed. This as- 
sumption was, relaxed in a recent paperEZI, where follow- 
ing our worka, the possibility of a generalized dynamic 
SDW in regions larger than the non-superconducting 
core, and co-existing with well-established superconduc- 
tivity, was appreciated. Hu and ZhangcJ also suggested 
that a small proximity-type coupling between the mag- 
netic domains centered on the neighboring vortices may 
be sufficient to stabilize static long-range magnetic or- 
der in a SC+SDW phase, in which enhancement of the 
SDW in the vortex cores was the dominant effect. As 
we have reviewed at length above, the strongly relevant 
exciton self-interactions lead to a different description of 
the SC+SDW phaseEl; in the SC phase, as one approaches 
the SC+SDW phase, the SDW order is induced over large 
length scales outside the vortex core, and the influence 
of the superflow is paramount. Only in the regime where 
magnetic field is small and the system is well within the 
SC phase (and far from the SC to SC+SDW boundary), 
can we speak in terms of localized bound state pulled 
below the continuum. However, this limit is of little ex- 
perimental interest, since it corresponds to high energy 
magnetic excitons (which, as discussed above, are proba- 
bly strongly overdamped by other excitations associated 



with the vortex cores) with a vanishingly small intensity. 



1. Pinning of charge order in the SC phase 

Our physical discussion has so far neglec ted t he in- 
fluence of the pinning potential in <Si at in 



1.12|) . We 



will continue to neglect this term in most of this paper, 
apart from computations in Section IV E whose content 
we briefly describe here. This analysis-is motivated by 
the STM experiments of Hoffman et al&B. 

The SC phase of Fig | preserves spin rotation invari- 
ance, and so has (<fr xa ) = and, by ( |1.2| ), (S a ) = (if we 
were to account for the small Zeeman term (Appendix^), 
the analogous statement holds for the spin density in the 
plane perpendicular to the magnetic field). In the ab- 
sence of 5i a t, all the remaining terms in the partition 
function Z in (1. 13) are invariant under the sliding sym- 



met ry <& xa — > e $ xol , and so we also have (& xa ) = and, 
by (1.5,Q|), (8p) — in the SC phase. Now if we include 
the effect of <Si a t perturbatively (which is all we shall do 
here), the pinning of the dynamic fluctuations by the 
vortex cores leads to static CDW order with (<& xa ) ^ 
and (dp) ^ 0, while the continued preservation of spin 
rotation invariance implies that we still have ($ 1Q ) = 
and (S a ) = 0. (Of course in the other SC+SDW phase, 
spin rotation symmetry is broken, and so (f m ) ^ and 
(S a ) + 1 0, along with static CDW order.) 

The nucleation of static CDW order, but with dynamic 
SDW order, in the SC phase by the vortices was first pre- 
dicted in Refs. [l0|5f], where a connection was also made 
with lattice scale studies of bond-centered charge order 
correlations in superconductors with preserved spin rota- 
tion invarianceES. These latter works found a significant 
doping range over which the charge order had a period 
pinned at four lattice spacings, which is the period ob- 
served in the STM experiments of Hoffman et al£3 (the 
same period also appeared in density matrix panormal- 
ization group studies by White and ScalapincEl) . Here 
we are interested in the spatial extent of the envelope of 
the period four charge order. Following Ref. [H], here we 
will compute this envelope using our present models for 
dynamic SDW/CDW fluctuations in the SC phase, and 
the pinning of a static CDW by <Si at . 

After this paper was originally released, we learnt of 
the microscopic model of Chen and TingH for the STM 
experiments, which follows the earlier work of Ref. fig . 
Their model has static order for both the SDW and CDW, 
and thus would apply only in the SC+SDW phase of our 
phase diagram in Fig [|. ft appears unlikely to us that 
the slightly overdoped BSCCO sample used by Hoffman 
et alH is in the SC+SDW phase. 

The simple model of the field-induced dynamic SDW 
fluctuations we have described in this section can be read- 
ily extended to compute the static CDW order induced 
by £ in the SC phase. Indeed, the upshot of our pro- 
ceeding discussion of the extended structure of So(r) is 
that we can use the Gaussian theory Sg in (1.15) with 
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V(r) given by the constant value in (1.21): computing 
($^„) in the theory Sq +&at f° r this value of V(r) and 
to first-order in (, we findtll 



<*2a(r,r)> = £ 



87r 3 / 2 [s eff (i?)] 1 /4 u 5/i 



-2\r-r„\y/s c!f (H)/v 
X ; ^ ,(1.23) 



V 



3/2 



,\= v ({{x - x v )/v x ) 2 + {(y - y v )/v 2 f) 



1/2 



where |r 

and v — {v\V2) 1 ^ 2 ] the result (1.23) holds for large |r— r v \, 
and the divergence at small |r — r v \ it cutoff by lattice 
scale effects. Note that the static CDW order decays ex- 
ponentially around each vortex core over a length scale 
= w/(2-v/s e jf (H)) which has been increased by the in- 
fluence of the field-induced superflow (by the decrease 
of s c g(H) in ( 1.20| )). Note also that this length scale is 
not related to any localization scale associated with the 
SDW state So(r); indeed, we have argued above that the 
latter state is extended. In the present simple Gaussian 
calculatio n, w e used the very simple constant potential 
given in (L21) in the Schrodinger equation for the ex- 



citon, (1.17); all eigenstates of such an equation are ex- 
ten ded p lane-wave states. Instead, the exponential decay 



1.23) arises from the integral over all the oscillating 



(but extended) excited states of (1.17). The body of the 
paper will show that the same feature also holds when 
the full form of V(r) is used, and not just the crude ap- 
proximation in ( 1.21 ) (see Figs |l5| and |l6|) . 

It is useful to make an analogy between the above re- 
sult and the phenomenon of Friedel oscillations in a Fermi 
liquid. A Fermi liquid state has no static SDW or CDW 
order, but there are enhanced fluctuations of these or- 
ders at 2kp, the wavevector which spans extremal points 
of the Fermi surface. In the presence of an external im- 
purity, static CDW oscillations at 2kp are induced, while 
full spin-rotation invariance is preserved. The amplitude 
of these oscillations decay with a power-law because the 
Fermi liquid has gapless spectrum of SDW/CDW excita- 
tions. 

In the present situation, the physics of the doped Mott 
insulator induces a preference for excitonic SDW fluctua- 
tions at the wavevectors K SX)JI and for CDW fluctuations 
at the wavevectors JCcx,y — 2K sa , )3/ . The SC phase has a 
spin gap, A, at these wavevectors, and so such spin corre- 
lations decay exponentially on the scale £ s = v/A (as we 
have noted, this is not a localization scale of the spin exci- 
ton states, which are all extended). The vortex core pins 
the phase the dynamic SDW fluctuations which reside 
above this spin-gap, and the resulting "Friedel oscilla- 
tions of the spin-gap" are manifested by static CDW os- 
cillations at the wavevectors K. cx ,y whose envelope decays 
exponentially over a length scale £ c = £ s /2. These may 
therefore be viewed as the Friedel oscillations associated 
with the excitonic bound states that are present below 
the particle-hole continuum of the d-wave superconduc- 
tor. In a weak-coupling BCS/RPA theory one can also 



expect additional Friedel oscillations associated with the 
continuum of particle-hole excitations, whose wavevector 
is determined by the quasiparticle dispersion. Such a pic- 
ture may be appropriate in the strongly overdoped limit 
in zero magnetic field, with pinning provided by impuri- 
ties. However, as one lowers the doping in the SC phase 
(to approach the boundary to the SC+SDW phase), an 
excitonic bound state appears, and we have focused on 
its physics here; the wavevector of this exciton is deter- 
mined by strong-coupling effects in the doped Mott insu- 
lator. The strength of this exciton could also be enhanced 
relative to the particle-hole continuum in the vicinity of 
vortices in an applied magnetic field-this effect requires 
explicit consideration of the fermionic quasiparticles, and 
so is beyond the scope of the theories considered here. 

The outline of the remainder of this paper is as follows. 
We will begin in Section || by a discussion of the phase 
diagram of the spin and charge density wave order param- 
eters in zero magnetic field. More complex phases and 
phase diagrams are also possible, associated the compos- 
ites and 'fractions' of these order parameters, but we will 
postpone their discussion until Section VI. We will turn 



to the influence of the magnetic field in Section III : here 



we will restrict our attention to the quantum transition 
described by <5>$, but most of the zero-field transitions 



discussed in Section VI have a related response to an ap- 
plied magnetic field. Section III contains a description of 



the phase diagram in the magnetic field, while the subse- 
quent sections describe the dynamic and static properties 
of the two phases on either side of the critical point in 
some detail: Section IV describes the SC phase, w hile 
Section [v] describes the SC+SDW phase. Section [VTI 



reviews earlier theoretical and experimental work on the 
interplay of magnetism and superconductivity, and dis- 
cuss its relationship to our treatment here. We conclude 
in Section VIII by considering implications of our results 
for recent experiments; readers not interested in theo- 
retical details may skip ahead to Section VIII now. A 



number of technical and numerical details appear in the 
appendices. 



II. PHASE DIAGRAM IN ZERO MAGNETIC 
FIELD 



We orient ourselves by discussing the phase diagrams 
of models with various types of spin and charge density 
wave order. We will restrict our attention in this section 
to zero external field, assume that a background SC order 
is always present in all the phases. As we have argued 
above, this implies that we need not consider the SC 
order parameter explicitly, and its influence only serves 
to renormalize various couplings in the effective actions. 
A somewhat different viewpoint, with a more explicit role 
for the SC order, has been taken recently by LeetJ. 

Here, we consider phases that are characterized simply 
by the condensates of one or more of the order parame- 
ters <& x ya and <j) x y , introduced in Section [A. More com- 
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FIG. 4: Mean field, zero temperature phase diagram of 
the zer o m agnetic field model 5$ + <S^ + 5$^ defined in 
( piljO]!^ , with u 2 < 0. 



plex phases associated with composites or "fractions" of 
these fields are also possible and these will be considered 
later in Section VI. However, the remainder of the paper 



will only deal with the influence of the magnetic field on 
phases and phase boundaries associated with the order 
parameters $ liS „ and 4> x ,y] the more complex cases have 
similar properties which can be described in an analogous 
manner. 

To characterize the simple phases we need an effe ctiv e 
action Sa, for the 



t> x ,yt while that for <& x ,ya is S$ in (1 



the former c an b e written down using a reasoning similar 
to that for (1.8), and we obtain 



S 



d 2 rdr 



\d T <j) x \ +v 1 \d x 4> x \ +v 2 \d y (j> x 

|5 T ^| Z +Vf \d X (j) y \ 2 + v\ \dy(j)y\ 2 + S(\(I) X ^ + \(/)y\ 2 ) 

(2.1) 



+ y(l^| 4 + l^| 4 ) + ^i|^| 2 |^| 2 



The correspondence ( |l.6| ) implies that for K. cx = 2K. SX 
and K CJ/ = 2K sy the SDW and CDW order parameters 
are coupled by 



S<s>d> = —A / d rdr 



/>*<£> 2 

X XOc 



b y & ya + c.c] ; (2.2) 



without loss of generality, we can assume that the cou- 
pling A > 0. At the mean-field level, the properties of the 
quantum model S$ + Sa, + S$0 are essentially identical to 
the classical models considered by Zachar et al.c3 for spin 
and charge ordering transitions at non-zero temperature; 
so we can directly borrow their results, and a character- 
istic mean-field phase diagram is shown in Fig H. 



Next, we discuss the critical properties of the various 
second-order quantum transitions in Fig [|. 

Near the transition between phase II (SC+SDW) and 
the symmetric phase I (SC), the primary order param- 
eters are <fe XtVa - We can integrate out the non-critical 
(j) x y fields and this merely renormalizes the couplings in 
S$. So the theory S$ is the critical theory for this tran- 
sition at H = 0. This is a model of some complexity, and 
the universal critical properties of related simpler jbc 
are the focus of some debate in the literatureE3£3'L3'E 
these earlier results are briefly reviewed in Appendix 
These previous studies correspond to the case where $ m 
and are decoupled (wi — W2 — W3 = 0), and weakly 
first-order transitions are obtained in some cases. We will 
address the generalization of these previous analyses to 
the case of non-zero u>i,2,3 in future work. Here, we will 
be satisfied by considering the simplest, and most sym- 
metric, case of a second-order transition: for the special 
values v\ — V2, ui = wi, U2 = W2 = W3 = the model 
S$ has a 0(12) symmetry, and its properties are identi- 
cal to that of the (N = 12)-component <^ 4 theory S v to 
be described shortly below. The influence of H on other 
second-order or weakly first-order transitions should be 
very similar, with the changes only modifying the nu- 
merical values of certain asymptotic critical parameters. 
Part of our reason for not expending much effort on this 
point is that these asymptotic critical are not particu- 
larly relevant for the experimental situation in H ^ 
anyway: after including the small effects of <Si a t in (1.12), 
the "sliding" symmetry of 5$ disappears, and the asymp- 
totic critical properties of the SC+SDW to SC transition 
in H 7^ become identical to the (N = 3) component 
<p 4 theory S v . We will discuss the H ^ properties 
of S v at some length in this paper, and we expect that 
closely related results apply to the generalized 5$ and to 
5$ + <Si a t • 

Near the transition between phases III and I in Fig ||, 
the roles of <& x ,yu and 4> Xty are reversed. Now we can 
integrate out the non-critical & x . yai this renormalizes the 
couplings in S^, and the renormalized S$ is the critical 
theory for this transition at H = 0. At non-zero H, a 
model closely related to the one discussed above applies. 
We will not explicitly present the results for this model 
here, as most physical properties are essentially identical 
to those of S$ + S^ + S$^ . 

The remaining second order quantum transition in 
Fig ^ is that between phases II and III. Both these phases 
have (4> x ,y) 7^ 0, and the charge order can be viewed 
as a non-critical spectator to the transition. For speci- 
ficity, let us assume that {<p x ) is real and positive, while 
(4> y ) — 0; other cases lead to similar final results. Now 
replace <p x . v by their expectation values in S$ + S$0 in 
(M), and examine fluctuations of $> XtVa at the Gaus- 
sian level: those of Re [$ xa ] have an energy lower than 
all other components. Close to phase boundary between 
II and III we can therefore assume that the critical the- 
ory involves only ^„(r, r) = Re [$ aa (r, t)], and all other 
components only renormalize the couplings in its effective 
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action. In this manner, we can conclude that the II to 
III phase transition is described by the familiar (N = 3)- 
component ip 4 field theory, with effective action 



5 K 



ePrdrj- [(d r cp a f + v 2 {V r <p a ) 2 

+ { S + n\^{v)\ 2 ) V l] + ^) 2 }, (2.3) 



where the index a = 1...N, and the field tp a (r,T) is 
real. We have rescaled spatial co-ordinates to make the 
velocities i»i,2 equal to the common value v. For com- 
pleteness, we have also includ ed th e coupling to the SC 
order ip which derives from (1 
prop erties of the theory J-/T- 



1.10) 



An analysis of the 
tp, defined in ( |1.9[ ) and 
(2.3), in non-zero field shall occupy us in most of the re- 
mainder of the paper. Recall also that the N = 12 case 
of this theory also describes a particular case of the I to 
II transition discussed earlier. 



III. PHASE DIAGRAM IN A MAGNETIC 
FIELD 



We now embark on a presentation of the main new re- 
sults of this paper: a description of the phase diagram 
and the dynamic spin spectra of T jT + 5„, defined in 



(1.9) and (2.3), as a function of the applied field H. As 
discussed near (2J3), this theory describes the response 
of a number of specific phase boundaries of states with 
SDW/CDW order to an applied magnetic field; the num- 
ber of components of (p a takes the values N — 3, 12 de- 
pending upon the transition of interest, but we expect 
similar results for all values of N > 3. Actually, closely 
related analyses can be applied to most of the phases 
to be discussed in Section [Vl| The basic effect, that all 
couplings associated with the non-superconducting order 
parameter acquire a Hhi(l/H) depends, is very robust 
and leads to analogous phase diagrams in almost all cases. 

The theory J-/T + has a rather number of cou- 
pling constants, and it is useful to use our freedom to 
rescale lengths, times, and field scales to obtain an irre- 
ducible set of parameters whose values control the struc- 
ture of our results. First, as is conventional in the stan- 
dard Ginzburg-Landau theory of superconductivity, we 
introduce the superconducting coherence length, £o, and 
the field scales H c and H® 2 '- 



H r = 



H° c2 




(3.1) 



as we noted earlier, H® 2 is the value of the upper crit- 
ical field at the point M in Fig ||, and H® 2 = V%kH c 
where n is usual the Ginzburg Landau parameter. Wo 



will also see below in Section IV B that the coupling a 
acquires a shift renormalization due to its coupling to (p a 
fluctuations: we assume that renor mali zation has already 
been performed in the definitions (3.1). We now use the 
length £07 the velocity v, and the parameters in (3.1) to 
set various length, time, temperature, field, and coupling 
constant scales; we define the dimensionless parameters 



H 





T = 




H 


if = 


ft 




V OL 




U = 


V' 


V 



T 



"0 ; fa = V V^oPa 



(3.2) 



It is evident from the above that we are measuring length 
scales in units of £o a- n d energy scales in units of v/£q. 

Collecting all the transformations, let us restate the 
problem we are going to solve; we drop all the tildes, and 
it is henceforth assumed that all parameters have been 
modified as in (p.2j). The partition function in (1.13) is 
now simplified to 



(3.3) 



where S v is as in (2.3) but with v — 1, while T is now 
given by 



T 



d 2 r 



-H 2 + ^| 4 + l(V r -*A)^| 2 



The dimensionless constant T is given by 

H 2 e d 



T = 



Attv 



(3.4) 



(3.5) 



where d is the inter-layer spacing (this factor of d is 
needed to make T dimensionless, and arises because F 
is the free energy per layer); in determining T, a useful 
unit of conversion is 1 (Jesla) 2 = 0.0624 meV A" 3 . The 
l in (|3.4|) now satisfies 



vector potential A 



V r x A = Hz. 



(3.6) 



An important property of the continuum theory ( [3.3] ) is 
that all dependence on the short distance cutoff can be re- 
moved by a single "mass renormalization" : this amounts 
to measuring the tuning parameter s in terms of its devia- 
tion from s = s c , the critical point between the SC+SDW 
and SC phases at H — 0. Consequently all physical prop- 
erties are functions only of the dimensionless parameters 
u, k, T, H/H® 2 , and s — s c . We will present numerical 
results for the frequency and spatial dependence of vari- 
ous observables below as a function of H/H® 2 and s — s c 
for the simple set of values u = K = T = 1; we do not 
expect any qualitative changes for other values of these 
last three parameters. Also, it will occasionally be con- 
venient to exchange the parameter s — s c for A, the value 
of the spin gap in the s > s c SC phase at H = 0. 
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The technical tool we shall use in our analysis of (3.2) 
is the large N expansion. This approachtJ is known to 
yield an accurate description of the vicinity of spin or- 
dering quantum critical points in two dimensions, and we 
expect the same to hold here in the presence of a non-zero 
H . Details of the approach will emerge in the following 
sections: here we summarize the main N = oo results for 
the positions of the phase boundaries appearing in Fig ^: 



• The tetra-critical point M where all four phases 
meet is at H = 1, s — s c = k. 

• The line BM represents the upper-critical held for 
the vanishing of superconductivity in the presence 
of SDW order; it is at 



A. Large N saddle point equations 



The index a in S v in ( |2.3| ) extends over a — 1 ... N , 
and depending upon the transition in Fig |] we are inter- 
ested in, we have either iV = 3orA = 12. For both cases, 
it is known that an accurate description of the physical 
properties is described by the large N expansion, whose 
implementation we shall now describe. 

First, we introduce an auxiliary held 

V(r,r) = S + K \^ H (r)\ 2 + 2u l pl(r 1 r). (4.1) 

We will often place a subscript H on various quantities 
(as for ip above) to emphasize that they are being evalu- 
ated at a non-zero H . Let us also denote 



H = 1 



+ 77^( s - s c)- 



4uT 4wT 



(3.7) 



• The line CM, the boundary for SDW order in the 
insulator, is at s — s c = k. 

• The line DM, the upper-critical held for supercon- 
ductivity in the absence of SDW order is at 



H = 1 



Nu 



8ttT 



'N 2 u 2 \V2 Nu 



/N 2 u- 
V 1~6^2 



(3.8) 



• Experimentally, the most important and accessible 
phase boundary is AM, the line representing onset 
of SDW order in the SC phase. The position of 
this line cannot be determined analytically: we will 
present detailed numerical results and an expansion 
in the vicinity of M; for small H its location behaves 
as 



H 



2(s - s c ) 
«ln(l/(s - s c )) ' 



(3.9) 



as may be readily dedu ced from (1.20), and was 
quoted already in ( 1.22j ). 



Our numerical as well as analytical studies will be di- 
vided into two parts, one for "SC" region of Fig || in 
Section and the other for "SC + SDW" region in 
Section 



IV. PHYSICAL PROPERTIES OF THE SC 
PHASE 

This section will describe an analysis of in the 

regime where spin rotation invariance is preserved with 
(ip a ) — 0. As we discussed earlier at the end of Sec- 
tion |l Cj, up on including the effect of the lattice pinning 
term (1.12) in a non-zero H, this phase does have static 



CDW order with (4> x ,y) 7^ 0, while preserving spin rota- 



tion invariance: this will be discussed in Section IV E 



s + K|V ff (r)| s 



(4.2) 



Now we add an innocuous term to S v , whose only effect 
is to multiply the partition function by a constant after 
a functional integration over V(r, r): 



r r 1/T 1 

S v -> S v - j d 2 r J dr — (V-2utp, 



,1/T 

d 2 r / dr 
Jo 

1 



-(d T <p a ) 2 + i(V r (^ Q ) 2 



V z 

8u 



(4.3) 



After integrating out ip a (a = 1 • • • N) , we have 

N, 



PV(r)exp 



-Trln(-a£- V; + V) 



-^Vs' + ^V 2 
4u 8u 



(4.4) 



Now by taking N — * 00 while keeping Nu constant, we 
obtain the saddle point equation in which V is a function 
of r but independent of r: 

V ff (r) = s + K\^ H (r)\ 2 + 2NuTY / G H {r,r 1 u; n ). (4.5) 

where the tp a propagator G#(r, r', u n ) is given by 

G H (v, r',cj n ) = (r| (a£ - V 2 T + V H (r))^ |r'), (4.6) 

with uj n a Matsubara frequency. In this case, the large- 
TV expansion is equivalent to a self-consistent one-loop 
calculation. 

The saddle point equat ion f or superconducting order 
parameter follows from ( 1.14 ): it is just the conven- 
tional Ginzburg-Landau equation with one additional 
term from the ip, ip coupling, 

1 + ^T^G ff (r,r,c„) 

+ |^(r)| 2 -(V r -zA) 2 l^(r) = 0. (4.7) 
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So the two unknown functions, Vf/(r), and ^g(r) are 
to b e de termined simultaneously by the solution of ( |4.5| ) 
and (4/7). As stated above, the expressions in these equa- 
tio ns dep end upon the short distance cutoff, but we show 
in ( [V B| ) that this can easily be removed by a simple shift 
of parameters. 



B. Renormalization of parameters 

It is first useful to obtain the complete solution of ( |4.5| ) 
and (LI) at H = 0. Let s = s c be the point where 
magnetic order appea rs (s o that (ip a ) ^ for s < s c ), 
whereB V = 0. Then (jij) tells us that 



/ dui f d 2 k 1 fJ . 
= s c + k\iJj Qc \ 2 + 2Nu — - (4.8) 



27T / 47T 2 W 2 +fc 2 ' 



where ipo c is the r ind ependent value of ip(r) at s = s c 
and H = 0, while (4.7) gives 



Nk f dw f d 2 k 



-i + liM' + ■==■ / -kz 



2T 



1 



2tt / 4tt 2 



0. (4.9) 



It is useful to normalize things so that ipo c = 1 at s = s c , 
H = and T = 0. This is achieved if we renormalize a 
to rem ove the offending term in (H.9|) . We make the shift 
in (l.S) (before the rescalings in fl3.2|)) 



a — > a 



/" da; /* d 2 /c 



2 a 



1 



27T / 47T 2 LU 2 + V 2 k 2 ' 



(4.10) 



Then, after (BJ), (4.7) is modified to 



f iV^r / /" dujd 2 k 1 

{- 1 + ^[ T E G *(^)-/^^ 



+ |^(r)| 2 ~(V r -zA) 2 }^(r) 



= 



k 2 

(4.11 



while (4.9) becomes simply 

V'oc = 1. (4.12) 
Now move to s > s c , where we have a spin gap 

(4.13) 



A = y/Vo > 0. 



Subtracting (4.8) from (4J5) we get 



An = s — s c 



■ k{\i/jo\ 



NuA Q 
2tt 



where (4.11) yields 



|^o| 2 = l + 



NkA 
8ttT ■ 



Inserting ( 4.15 ) back into ( 4.14 ) we obtain 



Nu 
~2k 



1 - 



K 2 \ 



4wT 



A = s - s c 



(4.14) 



(4.15) 



(4.16) 



Let us now use the above equations to simplify the 
equations for H ^ and T ^ 0. The new form will be 
independent of lattice cutoff. 

From (P~5jpl,pi~16l) we obtain 



V H (r)= A 2 + K (|V^(r)| 2 -|^ | 2 ) 
+2Nu[Tj2GH{r,r,uj n ) 

du d 2 k 1 



27T 47T 2 i 



fc 2 + A 2 



(4.17) 



where |^o| 2 is given in (4.15). Using (4.17) and (4.11 
obtain 



W'V 



1 - 



(|^(r)| 2 -|^| 



K 



(Vh(v)~A 2 ) 

{V r -iA)' z ^ ff (r) =0 (4.18) 



AuT J Vl " v 71 " U| ' 4uT 

i2 



The expressions (4.17 4. IS) are the main equations we 
shall solve for the unknowns V//(r) and ^ij(r) in this 
paper. It can be checked that at H = 0, T = 0, these 
equations are solved by Vjf = A 2 , and i/jh — ipo- We de- 
scribe the numerical solution of these equations for H ^ 
in Appendix [e] and present the results in the following 
subsection. A useful step in this numerical solution is 
the following pa ram eterization of the Green's function 
G fl (r,r>„) in (O) 



G H (r, r',w n ) 



E 



d 2 fc S; k (r)S Mk (rO 
1BZ 4^ 2 ^ 2 +£ 2 (k) ' 



(4.19) 



where ^^(r) are the complete set of c igcnf unctions of 
the analog of the Schrodinger equation (1.17) 



-V 2 + Vff(r))S Mk (r)=£ 2 (k)S Mk (r) 



(4.20) 



Here k is a 'Bloch' momentum which extends over the 
first Brillouin zone of the vortex lattice, /i is a 'band' 
index, and (k) are the corresponding energy eigenval- 
ues. All of our numerical analysis was performed for the 
values u = k = T = 1 and N = 3. 



Phase boundaries 



The equations (4.17) and (4.18) can be readily solved 
to obtain the locations of the CM and DM phase bound- 
aries in Fig H On DM, the superconducting phase pa- 
rameter ipjj (r) vanish es an d all parameters become r in- 
dependent; thus Eqn ( 4.17 ) becomes 



A 2 



Nu 

+ 2^ 
Nu 
"2V 



WV H ~ An) 



(4.21) 
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FIG. 5: The lowest eigenvalue of fll.l7| ), E (0) vs. H for 
s — s c = 0.2. The linear continuation of the line to solve 
Eq(0) — gives us the critical H for this s , which is about 
0.117 with an uncertainty of ±0.002. 




-0.5 
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FIG. 6: Numerical results for the phase boundary AM in 
Fig ^ for u = k = T = 1. Also shown is a portion of the phase 
boundary BM whose position is known analytically from (3.7). 



Different aspects of the physical properties are described in 
the remainder of the paper at the points labeled a-l. 



where we used Eqn ( [4. 15| ) and ( 4.14 ). Then from Eqn. 
(4.18) we have 



1 + 





K 2 


K 


1 - 


4uT + 


4uT 






K 


1 - 


4uT + 


4uT 


Nk 


^N 2 u 2 


— K 


8vrT 


A 167T 2 



(s -s c - V H ) 

< Nu /77- 
Z7T 

1/2 



Nu 

47T 



, (4.22) 



which is the result quoted in (3.8). Similarly, it is easy 
to see that the phase boundary CM is at s — s c = k. 

It remains to determine the location of the phase 
boundary AM, which is also physically the most inter- 
esting one. We dete rmin ed thi s bou ndary by a full nu- 



merical solution of ( 4.17 ) and ( 4.18 ) for a range of pa- 
rameters. Stability of t he S C phase requires that all the 
eigenvalues, E 2 (k), of ( 1.17 ) remain positive. The lowest 
of these eigenvalues is -Eb(O) and we followed its behavior 
as a function H: a typical result is shown in Fig [|. We 
expect Eq(0) to vanish linearly in the deviation from the 
critical field, as the critical theory is expected to be in 
the universality class of the ordinary 0(3) ip 4 field theory, 
and the latter has critical exponent zv = 1 in the large 
N limit. So we can determine the critical field by a linear 
extrapolation, and this is also shown in Fig ||. Combin- 
ing the results of such calculations at a range of values of 
s, we obtain our numerical result for the location of the 
AM boundary shown in Fig ^. 

Some further analytic results on the location of the 
AM phase boundary can be obtained in the vicinity of 
the multi-critical point M. It can be shown that the devi- 
ation of the phase boundary from M is linear in the large 
N limit i.e. it is at H = 1 — q(k — s + s c ), where g is 
a numerical constant. We describe these results in Ap- 
pendix |e], including the determination of g. The results 



obtained in this manner are consistent with our complete 
numerical analysis described above, and this is a strong 
check on our numerical analysis. 

Finally, we recall our result ( |3.9[ ) for the behavior of 
AM at small H and s — s c . Here there is a crucial log- 
arithm which follows from ( 1.20| ), and whose physical 
origin was discussed in Section [C. The signal of this 
logarithm are clearly visible in the phase boundary in 
Fig | 



D. Dynamic spin susceptibility 

In this section we describe the evolution of the dynamic 
spin fluctuation spectrum in the SC phase of Fig pi This 
is c lear ly specified by the Green's function Gjy(r, r', oj n ) 
in (4.6) which we computed above in determi nin g the 
phase boundary. More specifically, we see from ( tL.Sj ) that 
the observed dynamic spin susceptibility x(q, oj) is given 
by 



x(q, w) oc x v (q + K S2; , u) + x v (q - K 
+ X v (q + K sa , uj) + x v (q - K sy , 



sx, W) 

w), (4.23) 



where Xv> the dynamic susceptibility for the field tp a , is 
given by 



,(k,uO = i J d 2 rdVe ik -( r - r ')<3 H (r,r',( 
V / d 2 p S(p + G - k) 

,. /- JlBZ 



|c m g(p)| 
E 2 (P)-< 



(4.24) 



where V is the volume of the system, the p integration 
is over the first Brillouin zone of the reciprocal vortex 
lattice, G extends over the reciprocal lattice vectors of 



15 




point f 



FIG. 7: Reciprocal lattice of the vortex lattice. The density 
plots in Figs p[ [| [l^, and [n] have k varying along the arrow 
shown, with numerical values as shown. 



the vortex lattice, E^(p) are the eigenvalues of ( 4.20| ) (see 
als o Ap pendix [e]) , and the parameters c m g (p) ar e defined 
in (|E2|). We present resul ts for Im[x^(k, lo)] below. 

It is clear from ( 4.24 ) that in the present large N 
approximation, the spectrum of x<p consists entirely of 
sharp delta functions. These specify the dispersion of 
5 = 1 "excitons" which describe the SDW fluctuations, 
and are connected with the zero field "resonance" peak 
discussed early on in Section |. The excitons scatter off 
the vortex lattice, and our results describe the evolution 
of the resulting spectrum as one moves towards the onset 
of SDW order by increasing the applied magnetic field. 
We show the structure of Im[x v (k, oj)] by broadening the 
delta functions into sharp Lorentzians, and displaying the 
results in density plots. The momentum k in these plots 
varies along the direction of the reciprocal lattice shown 
in Fig[7. The results for a smaller value of s—s c are shown 
in Figs_| and 0, and those for larger value of s — s c are in 
Figs [10, and [Ty. Note that for very small H, there is 
less dispersion for the lowest mode: this is an indication 
that this excitation is centered on the vortex core, and 
there is weaker coupling between neighboring vortices. 
As the field is increased, this coupling increases, and the 
dispersion looks closer to that of a nearly-free particle, 
with weak reflections at the Brillouin zone boundaries of 
the vortex lattice. Also, the energy of the minimum ex- 
citation decreases with increasing field, until it vanishes 
at the AM phase boundary to the SC+SDW phase. 

We also show in Fig [l2]the spatial structure of the mod- 
ulus of the superconducting order parameter ^//(r)! 2 . 

The Brillouin zone boundary reflections above arise 
from the scattering of the exciton off the potential cre- 
ated by |^(r)| 2 . 

Finally, for experimental comparisons, it is useful to 
plot the intensity of the lowest exciton mode as a func- 
tion of the applied field. From ( 4.24| ) we see that this 
intensity is |coo(0)| 2 . We show a plot of this quantity in 
Fig |l4]. Observe that except for very small values of H, 
the intensity is of order unity, which is the behavior ex- 
pected for an extended exciton scattering off a periodic 



ntensity 
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FIG. 8: Density plot of lmx v in (4.24) in the SC phase for 
momenta along the arrow in FigH. The plot is for s — s c = 0.1 
and H — 0.01 (point / in Fig m). In this, an d all subsequent 
plots of \m\ v , the delta function peaks in (4.24) have been 
broadened into Lorentzians with energy width 0.01 for display 
purposes only. 
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FIG. 9: As in Fig g| but for larger H = 0.035, which brings the 
system very close to the AM phase boundary to the SC+SDW 
phase (point e in Fig []). 



potential as in Fig [T^. As H — ► 0, the behavior crosses 
over to that expected when the vortex cores are essen- 
tially decoupled, and the lowest mode is associated with 
a state loca.lized around each vortex core: in this limit, 
we expectEH the intensity ~ H . 
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FIG. 10: As in Fig |8j but with larger s - s c : H = 0.1 and 
s — s c = 0.9 (point m in Fig h). 
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FIG. 11: As in Fig [u] but with a large H = 0.8, which 
brings the system very close to the AM phase boundary to 
the SC+SDW phase (point Z in Fig |). 



E. Pinning of charge order 

This section will cons ider the consequences of the pin- 
ning term S\ a t in ( 1.12 ). We argued at the end of Sec- 
tion 1 C that this term pins the charge ord er, a nd leads 

^ (recall O)) in the 



to a static CDW with 
SC phase, while preserving spin rotation-invariance with 
(&x,ya) = 0. We have recently proposedHi this as an ex- 
planation for the CDW observed arou»d the vort ex i n the 
STM measurements of Hoffman et alx3. Section I C, also 



gave an initial estimate (in (1.23)) of the spatial structure 




FIG. 12: Spatial dependence of the modulus of the supercon- 
ducting order parameter \tpH (r)| 2 plotted on the rectangular 
half unit cell of the vortex lattice indicated by Fig |l3. This 
result is for s — s c — 0.5, and H — 0.1 (point k in Fig 3). 




FIG. 13: Half unit cell of the triangular vortex lattice in real 
space. 
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FIG. 14: Intensity of the lowest exciton mode in the SC phase, 
coo(0)| 2 , as a function of H for two values of s — s c . 
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of this pinned CDW: here we will obtain a more precise 
result, using the full solution of the SDW fluctuations in 
t he p resence of the vortex lattice. Using the relationship 
(1.6) between the CDW and SDW orders in the vicinity 
of the SC to SC+SDW transition, we conclude that to 
first order in ( 



0.2 n 



k,„(r)) cxCe- ra ^(r) 



(4.25) 



with 



f](r)ET^^G^(r,r„, Wn ), (4.26) 



where extends over the vortex lattice sites; clearly Q(r) 
has the full periodicity of the vortex lattice. 



We used our numerical solution of (4.17) and (4.18) to 
compute the function fi(r), which is proportional to the 
amplitude of the static CDW induced by the vortex lat- 
tice in the spin gap phase. We show our results for fi(r) 
in Figs [l5] and [l6] for a representative set of values in the 
SC phase. Also shown in the same figures, for orienta- 
tion, is the form of Soo (r) , the lowe st en ergy eigenfunc- 
tion of the dynamic SD W equ ation ( |l . 1 7| ) which appears 
in the Green's function (4.19). For very small field, both 
fi(r) and Soo(r) are localized around the vortex centers, 
with the localization length of the former being about 
half that of the latter. However, for larger fields, the 
exciton wavefunction Soo(r) gets delocalized, while the 
CDW order remains localized. This localization arises 



from the summation over all the states in ( 4.19 ) a nd is in 
keeping with the discussion at the end of Section I C 1 . 



V. PHYSICAL PROPERTIES OF THE SC+SDW 
PHASE 

W e now turn to the analysis of the partition function 
(3.3) in the phase with broken spin rotation invariance 
and (tp a ) = 0. This phase is reached when the lowest 
5 = 1 exciton mode in Section IV, S o(r), reaches zero 
energy (Eq(0) = 0) and then condenses. The presence 
of the condensate leads to long-range SDW order. We 
will adapt our large N computation to include such a 
condensate in the following subsection, and then describe 
the spatial structure of the condensate and the dynamic 
spin excitations. 



A. Large-TV saddle point equations 



The analysis here is parallel to that in Section IV A| . 
We introduce the auxiliary field V(r, r) defi ned in (4.1) 
and write the action in the form similar to (4.3). How- 
ever, to account for the condensate, we have to select a 
particular orientation in spin space, and treat the corre- 
sponding spin component in a selective manner. So we 
write 
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FIG. 15: Plots of the function Q(r) (filled symbols) in fl4.26p 
representing the static CDW order pinned by the vortices, 
along with the lowest SDW eigenfunction Eoo(r) of the dy- 
namic spin fluctuations above the spin gap (open symbols), at 
s — s c = 0.2. The spatial co-ordinate x is along the line con- 
necting two nearest-neighbor vortices and its scale has been 
chosen so that the vortex lattice spacing is unity (see Fig fL3h . 
The field takes the values H = 0.02 (squares, point h in Fign) 
and H = 0.1 (triangles, point g in Fig H); the latter field is 
close to the AM phase boundary in Fig B[ Note that the spin 
exciton state at point g is well extended through the lattice, 
while the charge order remains localized around the vortices. 
For point h the localization length of the spin exciton state is 
about twice that of the charge order . Th ese results are con- 
sistent with the discussion in Section FClL As was also noted 
below ( |l.23| ), the continuum expression ( |l.26| ) actually has a 
divergence for r equal to any r„: our numerical computation 
uses a finite momentum cutoff A, and this rounds out the di- 
vergence at distances |r — r„| < A -1 ; we have verified this by 
numerical computations at different A. In the same units as 
those for x in the figure, we used A ~ 36 above. 



and integrate out only 7Ti,7T2, • • ■ , 70v-i to obtain 



Z = J VV(r,T)Vn{r,T) 



exp 



-i-V S '-^(9 T n) 2 -^(V r n) 2 -^Vn 2 ], (5.2) 



(ViVn, 7n,7r 2 , • • • ,7Tjv-i), 



(5.1) 



where s' was defined in (4.2). Now we take N — > 00 while 
keeping Nu fixed, and ignoring the difference between N 
and N — 1. This leads to saddle point equations for the 
time-independent field V//(r) an d th e SDW condensate 
rt^f(r); these equations replace (4.5), but contain addi- 
tional terms due to the spontaneous spin polarization: 

Vff(r) = s + K\ip H (r)\ 2 

+ 2NuTj2G H {r,r,uj n ) + 2Nun 2 H (r) (5.3) 
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FIG. 16: As in Fig [TJ but for s - s c = 0.5. The field takes 
the values H = 0.1 (squares, point k in Fig 0) and H = 0.4 
(triangles, point j in Fig []); the latter field is close to the AM 
phase boundary in Fig |fj. Now both points have extended spin 
exciton states (that at point j is essentially a plane wave), 
while the charge order is exponentially localized. 



Next, subtracting (4.8) from (5.3) while noticing that i/joc 
is already renormalized to unity, we have 



Vff(r) = s - s c + «(|iMr)| 2 - 1) + 2Nun 2 H (r) 

dujd 2 k 1 



+ 2Ar,[T£GH( Wn )-/^g 



k 2 



(5.7) 



From (5.7) and (5.6) we have 



'-^l(l^r)P-l) + ^CMr)- S + Sc ) 



(V r - iA) 



ip H (r) = 0. 



(5.8) 



The final set of equ atio ns for the prop erties of the 
SC+SDW phase are Q, fTJ) and Q; these are to 
be solved for the unknowns Vjj(r), ^ij(r) and ns(r). We 
describe the numerical solution in Appendix |g| . 



Phase boundaries 



and 



(-V r 2 + V ff (r))n ff (r) = 0, 



(5.4) 



where Gh is given by ( |4.6| ). Comparing (5.4) and ( |4.6| ) 
it is easy to see that the spectrum of Gh, as defined in 
(4.19) has one mode with £o(k) — > as k — > 0; this is, 
of course, the Goldstone spin wave mode associated with 
the spontaneous SDW condensate. 

The equation which determined the superconducting 
order parameter, ?/>ij(r), was (4.7), and this is now re- 
placed by 



{-1 + -^[Tj2GH(r,r,u; n ) + 2Nun 2 H (r) 

+ |^(r)| 2 -(V r -zA) 2 }^(r) = 0. (5.5) 

B. Renormalization of parameters 



Now we pr ocee d as i n Secti on RV B to remove all de- 
pendence of (5.3), (5.4), and ( |5.5| ) on the short-distance 
cutoff. 

First consider the case when T = 0, H — and s = s r , 



where (4.8) and ( |4.9| ) hold. Now after we shift parameter 
a in as in ( 4.10 ), the Ginzburg-Landau equation ( |5.5| ) is 
modified to 



{ - l + ^[T£G ff (r,r,^ 



du> 
2^ 



d 2 k 



1 



47T 2 UJ 2 + k 2 

(V r -iA) 2 }^(r) =0 



2Nun 2 H {r) 
\Mr)\ 2 



We have already determined the posit ions of several 
phase boundaries in Fig || in Section IV C , and it remains 
only to determine BM. First notice that at the transition 
into a non-superconducting phase, the order parameter 
■0(r) vanishes, and thu s Vir(r) and Cf/( r ) are spatially 
uniform. So from (5.4) we have Vh = 0. Plugging this 



into (5.6) we obt ain the position of the phase boundary 



BM specified in (3.7) 



D. SDW order parameter 

The presence of a static spin condensate implies that 
the dynamic spin susceptibility contains sharp Bragg 
peaks at zero frequency and at wavevectors separated 
from the SDW ordering wavevectors by the reciprocal lat=. 
tice vectors of the vortex lattice as suggested by Zhango 
and discussed by us in Ref. |^; these are in addition to the 
dynamic spectra specified in ( f4.24 ). This means that the 
dynamic structure factor ^(k, uj) (which is related to the 



susceptibility x v (k, uj) in (4.24) by the usual fluctuation- 
dissipation theorem) has the contributions 

S v (k,w) = (2tt)*(w) \f G \ 2 (2TT) 2 S(k - G) (5.9) 

G 

where G extends over the reciprocal lattice vectors of the 
vortex lattice, and 



Au Ju 



d 2 i 



-iG-i 



(5.10) 



where the spatial integral is over U the unit cell o f the 
vortex lattice with area Ay. Note t hat, by ( 4.23Q , the 
physical momentum is related to k in (5.9) by shifts from 



(5.6) 

the SDW ordering wavevectors K. sx and K s „ 
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FIG. 17: Bragg scattering intensity, |/g| , as a function of 
H at s — s c = 0. Shown are the values at G = (squares) 
and at G = Gi (triangles), which is the smallest non-zero 
reciprocal lattice vector of the vortex lattice. Note that the 
intensities at G = Gi have been magnified by a factor of 150 
to make them visible on this plot. The intensities are zero at 
H — 0, because s — s c is the quantum critical point in zero 
field. The line shows 0.63 H ln( 3.61/_ff), which is the best fit 
to the functional form in (5.11). 




FIG. 18: As in Fig [jj but for s - s c = -0.3, showing |/ G | 2 
at G = (squares) and at G = Gi (diamonds) . Unlike 
Fig |l^, the intensity |/o| 2 is non-zero even at zero field. The 
intensities at G = Gi have now been magnified by a factor 
of 200. The line is 0.3 + 0.98iOn(2.12/#), which is the best 
fit to the functional form in ( 5.11 ) 



Figs |l7| and [T^ show plots of the Bragg scattering in- 
tensity |/g 1 2 , for the two smallest values of G and two 
values of s — s c , as a function of H. 

As argued in Ref. [| the correspondence (1.20) implies 
that the scattering intensity at zero wavevector, |/ | 2 
should increase with as 



(|/o| 2 ) cxHln(l/H). 



(5.11) 



The fits to this functional form in Fig [T?] show that this 
works quite well. Notice also that the intensity at the first 
non-zero reciprocal lattice vector, Gi, is quite small, and 




FIG. 19: Spatial form of the SDW order parameter nn(r) 
in the SC+SDW phase at s - s c = -0.3, H = 0.35 (point 
b in Fig |^) over vortex lattice shown in Fig ^| Notice that 
the vertical scale extends over a rather short range, and the 
modulation in ran(r) is quite small relative to the uniform 
component. 



that it decreases with increasing H. This suggests that 
observation of this satellite peak is best performed at 
as small a field as possible-of course, H should be large 
enough so that |Gi| is large enough to be outside the 
resolution window of the peak at G = 0. It is interesting 
to observe here that we can view the Bragg peak at Gi 
as arising from condensation at the non-zero k minimum 
in Fig [| of the dispersion of the exciton in the SC phase. 

For completeness, we also show the real space form of 
the condensate rifl-(r) in Figs [l9| and ^ for two points in 
the SC+SDW phase. The spatial form of the modulus 
of the superconducting order parameter for the first set 
of parameters is shown in Fig 
analog of Fig [ll| which was for the SC phase. 



21, This last figure is the 



E. Dynamic spin susceptibility 



Finally, we follow the presentation in Section IV D 
and discuss the dynamic spin spectrum in the SC+SDW 
phase. The non-zero u spectral densities prese nted here 
appear along with the u> = contributions in (5.S). We 
will restrict our attention to the susceptibility transverse 
to the ordering direction: this is given by the fluctua- 
tions of the last N — 1 components in (5.1), whic h are 
in turn related to the Green's function Gh m (5-3) and 
( |5.5p . So the transverse dynamic spin susceptibility is 
given by a formula analogous to (4.24). As before, we 
present the results by broadening the delta functions to 
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FIG. 20: As in Fig [t|; in the SC+SDW phase at s - s c = 0.5, 
H — 0.45 (point i in Fig []) over a single vortex lattice unit 
cell shown in Fig Fil 



point d 




FIG. 22: As in Fig g, but for the transverse susceptibility in 
the SC+SDW phase. The parameter values are s = s c and 
H = 0.05 (point d in Fig |). 



point c 




FIG. 21: Spatial dependence of the modulus of the supercon- 
ducting order parameter \-^ H {r)\ 2 plotted on the rectangular FIG. 23: As in Fig |5| but for the values a = s c and H = 0.3 
half unit cell of the vortex lattice indicated by Fig As in (point c in Fig []). 
Fig this result is for s = s c , and H = 0.35 (point b in 



sharp Lorentzians. 

Our results for the spectral densities are shown in 
Figs ^2| H^, and |25| for a series of values of s — s c 
and H in the SC+SDW phase. Note first that there 
is always a gapless spin-wave mode. In addition there 
are features arising from scattering off the vortex lattice: 
these are strongest in the vicinity of the quantum critical 
point A at zero field. 



VI. OTHER PHASES IN ZERO MAGNETIC 
FIELD 

The next two sections involve a slight detour from the 
main flow of the paper. This section we will examine 
phases and phase transitions associated with composites 
or 'fractions' of the primary order parameters $> x ,ya and 
4> Xt y. This is done mainly for completeness. Readers not 
interested in this detour may skip ahead to next section 
without loss of continuity. 
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FIG. 24: As in Fig but for the values s — s c 
H = 0.35 (point a in Fig |). 
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FIG. 25: As in Fig |22j but for the values s — s c = 0.5 and 
H = 0.45 (point i in Figg). 



Phases with nematic order 



In Section I_A we argued that a generalized non-two- 
sublattice spin density wave order may be associated with 
a charge density wave. Another interesting possibility is 
that of spin nematic order, which has been previously 
discussed in Refs. |6(||6l"l|62||63| . If the CDW order param- 
eter may be understood as a spin zero combination of 
two S a operators (Sp(r,r) ~ S^(r, r)) , then the spin 
nematic order parameter Q a p{v, r) corresponds to their 



Qc/3(r, r) 



Sa (r,T)S p (r,T)- S -fs 2 a (r,T) 



(6.1) 



We pause briefly to also mention here an "Ising na- 
matic" order which has also been considered recentlyo. 
This order resides in real space associated with the lat- 
tice, and is distinct from the spin-space nematic order we 
are considering here. Order parameters with the Ising 
nematic order are I^qI 2 — l^al 2 and \(f> x \ 2 — \<t>y\ 2 i and 
these clearly measure a spontaneous choice between the 
x and y directions of the lattice. Our effective actions for 
&x,ya and (j) x ^ y are rich enough to also allow such orders. 

R eturning to our discussion of spin nematic order in 
( |6.l|) , we see that spin nematic order para meters that 
are consistent with the SDW order in (1.2) may be at 
wavevectors (0,0) and "K cx ,y 



? Q /?(r,T) = Q a/?(r,T) 

-Re[Q XQ/3 (r,r) e iK - r - 



yap(r, r)e 



(6.2) 



It is natural to call Qoa/3 a uniform spin nematic or- 
der parameter, and Q x , ya i3 a spin nematic density wave 
(SNDW). Both order parameters are symmetric (Qi a p — 
Qif3a)i but the uniform spin nematic Qoq/3 must be real, 
and the spin nematic density wave Q x ,ya/3 m ay be com- 
plex. The uniform spin nematic couples to the SDW 
order parameters & x . ya as 



(6.3) 



The spin nematic density wave Q xa p(r, r) couples to <& xa 
via 

Sq*,$ x = - A 2 J d 2 rdr [ Ql a p ( $ xa $ x/3 - i 5 a[i $ 2 5 ) 

+ c.c. ] (6.4) 

with a similar coupling between Qyapir, r) and & ya - 

The effective action for the spin nematic order param- 
ete rs m ay be written from the analysis of the symmetries 
of (6.2). The interplay of the spin nematic and spin den- 
sity wave orders may produce an extremely rich phase 
diagram. We will not attempt to explore its full richness, 
but restrict ourselves to the discussion of some simple il- 
lustrative examples. It is also worth pointing out that the 
appearance of the spin nematic order (either uniform or 
SNDW) does not give rise to the additional Bragg peak 
at zero energy, but produces a difference in the scattering 
cross sections for different neutron polarizations. 



1. Uniform spin nematic 

To write the effective action for the uniform spin ne- 
matic Qoa/3 we can give essentially the same arguments 
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as in deriving the Landau free energy for the classical 
nematics (see e.g. Ref]64|) 



Sqo - 



d rdr 



( d T Qo a f3 ) (d T Qop a 



+ Vq( VQqq/3 ) ( VQo/3 Q ) + A QoapQopa 



1 



' C*2 Qoa/3Qoi3~fQo~<8QoSc 



(6.5) 



By an appropriate spin rotation, the uniform spin ne- 
matic order parameter may always be brought into the 
diagonal form (this follows from the fact that it is a real 
and symmetric matrix) 



Qoq/3 



-Uq + v) 




(6.6) 



When (q) =/= but (r/) = we have a uniaxial spin ne- 
matic, and when both expectation values are finite we 
have a biaxial spin nematic. 

Let us start by considering the interplay of the uni- 
form spin nematic with the collinear SDW (for simplic- 
ity we only consider one of the SDW orders, say 
A schematic mean-field phase diagram at T — for 
5$ + Sq + 5q ,$ with BXi < and 1*2 < is shown on 
Fig. Thick lines correspond to the first-order tran- 
sitions, and thin lines correspond to the second order 
transitions. Phase I (SC) has no magnetic order of any 
kind; phase II (SC+SDW) has commensurate SDW or- 
der, which is accompanied by a uniaxial spin nematic 
order; phase III (SC+UN) has a uniaxial spin nematic 
order. For BX\ > and 1*2 < the phase diagram qual- 
itatively remains the same, however the phase II has a 
finite expectation value of both q and 77 in (6.6), so it has 
an SDW order accompanied by the biaxial spin nematic 
order. A schematic phase diagram in the case U2 > 
is shown on Fig. The phase II may now be a circu- 
lar spiral SDW (Ha), an elliptic spiral SDW (lib), and a 
collinear SDW (lie). 
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FIG. 26: Mean field zero temperature phase diagram of the 
model <S* + Sq + Sq ,* in zero magnetic field for the case 

U 2 < 0. 
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FIG. 27: Mean field zero temperature phase diagram of the 
model <S$ + Sq + Sq 0i <s> in zero magnetic field for the case 

U 2 > 0. 



2. Spin nematic density wave 

For the spin nematic density wave the third order 
terms are prohibited by symmetry: they carry oscillating 
factors e ±lKcx '!' r , and vanish after integrating over space 
in the long wavelength limit. Hence, 



= / d 2 rdr 



(drQt a p) (d T Qxf3a) 



+ T Cl ( Q X xa aQxfia ) 2 + T C 2 Q^nQx^Q^Qx&a 



I ~ 



C*3 Q xaftQ x8 r yQ x ~ffiQ x fi° 



(6.7) 



and there is a similar action Sq . 

The order parameter for the spin nematic density wave 
can be conveniently written using five complex numbers 
(see also Ref. |65[) 



Qxafl 



/ 4>xi , 

lpx3 



V 



x4 



Ipxl 
■0x5 



x4 

+ 1px2 4>x5 



(6.8) 



V3 / 



with normalization condition X^ a =i 5 iV'zal 2 = !• This 
representation makes obvious the connection between the 
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order parameter for the spin nematic density wave and 
condensates of spin 2 particles, for which Ciobanu et.alEl 
argued that there exist three distinct phases (not related 
to each other by spin rotations), depending on the pa- 
rameters Ci, C2, and C3. 

The phase diagrams of the spin nematic density wave 
order vs the SDW order is similar to the case of uniform 
spin nematic (Figs. ^6] and p7| ) with the main difference 
that the phase boundary between I and III is now second 
order. 



the Z2 gauge transformation 



a, 

for <7j = ±1. 

Models of the ki 
in various contexts 63 " 6 ' 



7T 



(6.11) 



ie been discussed earlier 
o8HS9U7cU7i| j t wag p i n t e d ou t ; for 



example, that another term allowed by symmetry is a 
Maxwell term for the lattice gauge field 



B. Exciton Fractionalization 

Before concluding the section on the phases in zero 
field we would like to point out another interesting pos- 
sibility for the system described by the generalizations of 
5$. Consider this model in the regime where the spiral 
fluctuations are strongly suppressed, so we need to con- 
sider the collinear SDW order only; this happens in ( |1.S| ) 
for U2 < and with |tta| large. For simplicity we restrict 
our discussion to a SD W at wavevector K SJ) , & ya . As 
discussed in the Section I A, the collinear SDW can be 



written in the form in (1.3), which we reproduce here for 
completeness: 



$j, Q (r,r) = e l 



«(r, r). 



(6.9) 



We also noted below fll.3| ) that such a separation of the 
physical order parameter $ ya into the phase 9 and the 
real vector n a has an implicit ambiguity as we can simul- 
taneously change the sign of both without altering <fr ya . 
Formally this means that, for incommensurate K SJ/ , the 
order parameter Q ya belongs to the space (S2 X S\)/Z2- 
For commensurate K sa = 2717// (pa), where p' , p are 
relatively prim e integers, higher order terms not con- 
tained in ( |l.8| ) (but mentioned below it) imply 9 prefers 
a discrete set of valuesEHlEJ and the space is restricted to 
(52 x Z p )/Z2- Also, if full SU(2) spin rotation symmetry 
is absent, and the spins have an easy-plane restriction, 
then the first 5*2 factor changes to Si. 

The Z2 quotient in the order parameter space can be 
explicitly implemented as an Ising gauge symmetry, and 
it puts important constraints on the effective low energy 
theory. The lattice model consistent with such symmetry 
has the form 



Si = \_. J s <Jijn ia n,j a + 
(ij) 



^J c a lj cos{9 l - 0j) (6.10) 



where i and j are sites on the space-imaginary time lat- 
tice, the sum over (ij) extends over nearest neighbor 
links of this lattice, J s and J c are couplings imposing 
the propagation of SDW and CDW order respectively, 
nia = n a (Yi,Ti), &i = 0(rj,T,), and a i3 = ±1 is an Ising 
gauge field that lives on the links of the lattice. One can 
easily see that the lattice action (6.10) is invariant under 



(6.12) 



where the sum on □ extends over the plaquettes of a 2+1 
dimensional lattice. Such a term may be generated by in- 
tegrating out the high energy degrees of freedom or may 
be present due to certainpEcuattating terms in the original 
microscopic HamiltonianoEijO. This term has a striking 
effect on the properties of the model ( |6.10] ): it gives rise 
a phase in which the exciton $ yQ , fractionalizes, and fluc- 
tuations of n a are separated from the fluctuations of 9. 
Loosely speaking, the SDW and the CDW fluctuations 
get decoupled. 

It is useful to discuss the consequence of the 
confinement-deconnnement in the symmetric phase in 
which global symmetries are preserved: the models of this 
paper are invariant under SU(2) spin rotations, and the 
sliding U(l) symmetry (for commensurate values of K stf , 
the U(l) symmetry is reduced to a discrete Z p 'clock' 
symmetry, but essentially unchanged considerations ap- 
ply neverthelesai3lid) . The immediate manifestation of 
the confinement-deconfinement transition in such a sym- 
metric phase is the change in the degeneracy of the lowest 
energy excitations. In the confining phase their degen- 
eracy is 6: this 6-fold degenerate excitation corresponds 
to the quanta of the exciton field <fr ya , which have 6 real 
components. In contrast, in the deconfming phase we 
have separate excitations with degeneracies of 3 and 2, 
corresponding to quanta of rii a and 6i respectively. This 
may be understood by noting that the unbroken sym- 
metry ground state of the model Si + S a is a singlet 
ground state of the 5*0(3) x SO{2) rotors, where in the 
confining phase the angular momenta of the two rotors 
(h,h) = (Lso(3), Lso(2)) are bound by the constraint 
li+h = even, but this constraint is not present in the de- 
confining phase. Hence, in the confining phase the lowest 
excitation has (li = ±1, 12 — 1), which gives the degener- 
acy of 6. In the deconfming phase we can have excitations 
(li = ±1, 12 = 0) and (li = 0,h = 1), and these have de- 
generacies 3 and 2 respectively. We point out that the 
exact degeneracy of h = ±1 states requires the absence 
of the Berry's phase for the the 5*0(2) rotor, and comes 
from t he in version symmetry of the system, as was noted 
below (1.8). It is not related to the possible particle- hole 
symmetry of the underlying microscopic model. 
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It is worth emphasizing that the exciton fractionaliza- 
tion discussed above has a very different physical inter- 
pretation from that of electron fractionalizaiipn discussed 
in 'RVB' theories of doped Mott insulatorscj: in the lat- 
ter there are elementary S = 1/2 spinons which do not 
appear in our fractionalized states above. Instead our 
exciton fractionalization is within the sector of spin and 
charge density waves, and the collective spin excitations 
only have integen.spin. 

Zaanen et alx3 have recently discussed fractionaliza- 
tion in a microscopic picture of spin and charge order in 
"fluctuating stripe" states: the physical content of their 
analysis is quite similar to that of our discussion above. 
However their proposed effect action does not include the 
CDW phase field 6*.;, and we believe this is essential for a 
complete description of stripe physics. 

We have implicitly assumed above that the exciton 
fractionalization transition occurs in a background of 
SC order. However, a similar transition is also possible 
within a Fermi liquid. We believe that such a quantum 
critical point is a promising candidate for describing the 
finite temperature crossovers in the normal state of the 
cupratcs. Ordinary SDW/CDW transitions in a Fermi 
liquidEj have the unsatisfactory (in our view) feature of 
flowing to a free field fixed point because they are in 
their upper-critical dimensions. In contrast, the exci- 
ton fractionalization transition may well remain strongly 
coupled even in the presence of Fermi surface. Corre- 
sponding speculations of fractionalization influencing fi- 
nite temperature quantum criticality were also made by 
Zaanen et alE3. Again, their and our proposals should 
be distinguished from those associated with electron frac- 
tionalization made in e.q. Ref. |73- 



C. Topological defects 

An alternative picture of fractionalization, and of the 
various order parameters above, may be given in the lan- 
guage of the topological defects of the SDW phase; the 
condensation of distinct defects in the SDW state dis- 
tinguishes the new phases that appear. To simplify the 
presentation of this subsection we will describe the case 
of an easy plane antiferromagnet, in which the vectors 
&ia = e l6i ni a may only be in the x-y plane, but will also 
state the results for systems with full SU(2) spin rotation 
symmetry. A related discussion of defects in SDW states 
also appears in Ref. |7^. 

We start by giving a simple cartoonE3E3E3 of the non- 
two-sublattice SDW order & ya = const and the associ- 
ated CDW in Fig. Hole rich stripes (indicated by the 
dashed lines) act as antiphase domain walls for the hole 
poor antiferromagnetic domains. The Neel order shown 
by arrows changes sign when crossing such domain walls 
(the Neel order should not be confused with the vector 
n,i a which appears in the definition <F; Q = e l6i ni a ; the 
former oscillates as shown in Fig [2^, while rii a is con- 
stant in this configuration.). 



FIG. 28: A schematic picture of the non-two-sublattice 
collinear SDW order and associated CDW as a periodic array 
of anti-phase domain walls in Neel order at (it /a,n /a). Ar- 
rows show the change of sign of the Neel order across a hole 
rich domain wall. The fields m a and 9i are space-independent 
in the above configuration. 



Schematic pictures of the topological defects of the 
collinear SDW state are shown on Figures ( p9| ) - ( |3l| ) 
with crosses indicating the locations of the centers of de- 
fects (see also Ref. p2J) . These defects can also be for- 
mally classified by computing the homotopy groups; for 
systems with anpfiasy-plane spin symmetry the relevant 
homotopy grouplj is iri((Si x Si)/Z 2 ) = Z x Z, while 
for full SU(2) spin symmetry it is tti((5 2 x Si)/Z 2 ) = Z. 
These mathematical statements actually obscure some of 
the physical content, as will become clear from our dis- 
cussion below. 

We first discuss the physical content of the defect clas- 
sification for the easy-plane case. Consider the most ele- 
mentary topological defect: this is a composite of a 1/2 
vortex for the phase Oi and a 7r-disclination (i.e. 1/2 a 
meron) for the vector rii a (see Fig ^9|) ; this defecLis also a 
central actor in the discussion of Zaanen et alE3. When 
circling around such a defect both e lS and n; Q change 
sign, however the physical order parameter — e i ni a 
is uniquely defined. Given the circulations in 9i and 7ij ja , 
we label this defect (1/2,1/2). Actually, we can make 
four such elementary defects by changing the signs of the 
circulation of 9i and 7r-disclination and taking all of such 
combinations: we label these as (±1/2, ±1/2) in an ob- 
vious manner. Pairs of such elementary defects may be 
combined to give a full vortex for 8, which is trivial in 
the riia sector (see Fig 3^; this is the defect (1,0)) and 
a meron for the rii a , that it is trivial in the 6i sector 
(see Fig |3l]; this is the defect (0, 1)). Continuing in this 
manner, we see that all defects are labelled (mi/2, 7712/2) 
with mi , mi integers such that m,\ + to 2 is even. These 
labels lie on the analog of a FCC lattice in two dimen- 
sions; this is equivalent to a square lattice after a ro- 
tation by 45 degrees, and hence the homotopy group is 
Z x Z. This mathematical statement hides the fact that 
there is a fundamental physical difference between the 
(±1/2, ±1/2) and the (1,0), (0,1) defects, which we have 
discussed above. 

Next we turn to the case with full SU(2) symmetry. 
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FIG. 29: Elementary topological excitation of the collinear 
SDW phase: a composite of 1/2 vortex in 8i and 7r-disclination 
in riia . Both e l8 ' and n a change sign when going around this 
topological defect, but the physical order parameter <& ya i = 
e l6i rn a is single valued. 
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FIG. 30: Elementary topological excitation of the collinear 
SDW phase: a vortex in 9i. The circulation of 8 is equal to 
2tt. 
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FIG. 31: Elementary topological excitation of the collinear 
SDW phase: a meron of rii a . Such an object is stable only 
in systems with an easy-plane symmetry. Far away from the 
vortex core rii a winds in the plane like a usual vortex. Closer 
to the vortex center it may acquire an out of plane component. 
Systems with full SU(2) spin rotation symmetry only have 
point-like, instanton defects in spacetime: hedgehogs. 



dense, but 1/2 vortex - tt disclination composites 
remain gapful excitations, we find exciton fraction- 
alization as discussed above. The uncondensed 1/2 
vortex - 7r disclinatioijLcomposites correspond to the 
finite energy 'visons'Lj of the fractionalized phase 
of the Z2 gauge theory. 

• When only the merons (or hedgehogs) condense we 
find the CDW phase with no spin order. 

• When only the 9- vortices condense we get the spin 
nematic phase with no CDW order. 



VII. EARLIER WORK ON SC AND SDW 
ORDERS 



Now the 1/2-meron in rii a is actually equivalent to the 
— 1/2 meron (they are both better called 7r disclinations), 
and so there is no distinction between (1/2,1/2) and 
(1/2,-1/2); moreover, the (0,1) defect is topologically 
trivial. Consequently the spacetime line defects can sim- 
ply be labeled mi/2, where mi is an integer representing 
the phase winding of 9i, and hence the homotopy group 
is Z . However, there continues to be a fundamental phys- 
ical distinction between the cases with mi odd and even. 
For mi odd, there must be a corresponding tt disclina- 
tion in rii a , while for mi even the rii a configuration can 
be constant. The SU(2) case also has point defects in 
spacetime, the 'hedgehogs', which proliferate at spin dis- 
ordering transitions. 

The various phases discussed above can be easily un- 
derstood using the picture of topological defect conden- 
sation in a phase with conventional SDW order (the 
SC+SDW phase): 

• When the elementary 1/2 vortex - tt disclination 
composites condense we have a conventional (un- 
fractionalized) disordered phase (the SC phase). 

• When vortices and merons (or hedgehogs) con- 



As we noted earlier, this section is a detour from the 
main flow of arguments in this paper. For completeness, 
we review earlier theoretical and experimental work on 
the interplay of magnetism and superconductivity, and 
discuss connections to our treatment here. Less special- 
ized readers may skip ahead to the conclusions if they 

Early neutron scattering measurementsElIll of the evo- 
lution of the magnetic order in La2_a(Sr,EkL),5Cu04 with 
6 observed spectra which were interpretedo as evidence 
for the proximity of a quantum critical point at which the 
SDW order vanished, and which-aWwed strong hyperscal- 
ing properties. It was proposeccaO'E2l that such a quan- 
tum critical point (with dynamic exponent z = 1) con- 
trolled physical properties over a range of doping concen- 
trations. Further support for such a proposal appeared in 
the NMR experiments of Imai and collaboratorsEJ which 
displayed crossovers characteristic of the vicinity of a 
magnetic quantum critical point, with the critical point 
at a doping concentration 6 — 6 C ~ 0.12; similar evi- 
dence was presented recently by Fujiyama et alEn (for 
a review of the NMR data in this context, see Ref. [|). 
The concentration 6 C — 0.12 is well within the super- 
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conducting phase, and so the magnetic transition takes 
place within a background of superconducting order i.e. 
there is a second order transition between a phase with 
co-existing SC and SDW order (the SC+SDW phase) 
and an ordinary superconductor (the SC phase). As 
we noted in Section |lL, the neutron scattering measure- 
ments of Aeppli et ala at S — 0.14 provided rather di- 
rect evidence for such a magnetic quantum critical point. 
Additional evidence for microscopic co-existence of SC 
and SDW orders has appeared in a number of recent 
experimentsM'HHElfejBEaEZl. 

(For completeness, we also note here the additional 
phases present at very low 8 which were not the sub- 
ject of study in this paper: in La2_5Sr^Cu04, the 
three-dimensional, two-sublattice, insulating Neel state 
is present for 5 < 0.02, and is followed by an insulating 
SDW state with its wavevector polarized along the diag- 
onal (1,±1) directions™! As noted in Section Q, at 8 = 
0.055 there is a first-order insulator-to-superconductor 
transition to the SC+SDW phasefilfl, which has the SDW 
oriented along the (1,0), (0,1) directions; we discussed 
the properties of this SC+SDW phase in this paper.) 

A significant implication of the existence of a magnetic 
critical point at S = S c is that remnants of the magnetic 
excitations should be visible in the SC phase at 8 > 8 C . 
As originally discussed in Ref. Ej5|, for such critical points 
there is a sharp, gapped S — 1 collective mode (a spin 
exciton) which would appear as a 'resonance' in the neu- 
tron scattering cross-section. This resonance should ap- 
pear at the SDW ordering wavevector in (1.1), and recent 



evidence for gapped, low energy spin excitations at such 
a wavevector is in Refs. ^8|,^9[ Strong resonant scattering 
is alsi^een-al-the Neel order wavevector, (it, ir) in the SC 
phasdSS'L^LSEj: this remains at relatively high energies 
and may be viewed as a remnant of commensurate cor- 
relations at short length scaleal3. Batista et aLEHl have 
argued that the strong gapped response at (n, ir) is due to 
the superposition of the response at the two neighboring 
SDW ordering wavevectors at ±i? in (1.1). 

Another perspective on this quantum critical point, 
which was useful in our analysis, was provided by Zhang's 
SO (5) theoryo. This theory goes beyond the picture of 
competing SC and SDW orders in the ground state and 
adopts a stronger assumption of a microscopic dynamic 
symmet ry ,, bptsre cn them; this-Ja^s been supported by 
analyticr 4 lr 5 tr 6 H 97 l and numeric! 9 tt"! studies of a number 
of models. The generator of the enlarged SO (5) sym- 
metry is the 7r-excitation, a S =. _L. collective mode with 
charge 2 and momentum (ir, 7t)lj'L3. A sharp distinction 
between the models with and without the n excitation 
is possible in the weak interaction limit of a generalized 
BCS-RPA theory, where by going to the normal state one 
can check for the existence of a sharp collective mode 
with the quantum numbers of the 7r-particlal23. How- 
ever, a clear distinction is absent in the physically rel- 
evant strong coupling regime. For example, in the SC 
phase charge is only conserved modulo 2, and this charge 
2 particle is in fact indistinguishable from the neutral 



5 = 1 exciton in earlier theo ries^ of the SDW ordering 
transition (see also Ref 101). This exciton is smoothly 
connected to the S — 1 excitation in a paramagnetic 
Mott insulatorE23, and an interpretation of its 'resonance 
peak' as a generator of SO (5) rotations does not hold. In 
zero applied magnetie.iield, it is possible to formulate a 
theory of the excitonaEJ, and the associated SDW fluc- 
tuations, without any explicit reference to the SC order; 
the SC correlations only serve to modify various cou- 
plings in the effective action for the SD.W order. What 
we abstract from the analysis of Zhangc3 is the idea that 
the strength of the SC order itself should be viewed as 
a parameter which tunes the system across the magnetic 
quantum critical point: this emphasizes a local competi- 
tion between the SC and SDW orders. 

We also mention that these SO(5) models naturally 
describe a competition between the SC and the two sub- 
lattice SDW (Neel) phases. Non-two sublattice SDW can 
then appear as a result of the competition betwee n phas e 
separation and long range Coulomb interaction! 103 !! 104 !, 
across a first-order transition from the SC to the SDW 
phase. In this paper we will describe effective models for 
the non-two sublattice SDW directly, across a second- 
order transition from the SC to SC+SDW state. 

The precise nature of the interplay of SC and SDW 
orders in the cuprates at non-zero temperatures in three- 
dimensional models been a controversial subject (this 
paper dealt with two-dimensional quantum models at 
T = 0, and so the issues in this paragraph are only pe- 
ripherally related tcyxuir main discussion) . Following ear- 



lier general analyse d 105 ! Zhango pointed out four generic 
possibilities for the phase diagram, proposed the appear- 
ance of exact SO (5) symmetry in the classical theory of 
a finite temperature bi-critical point (this symmetey-|is 
actually only present in the equal-time correlatoralPJa) , 
and suggested that this is the situation most likely real- 
ized for the cuprates. In the presence of such a bi-critical 
point, there is a first-order transition between the SC 
and SC+SDW phases at low temperatures, and the en- 
ergy of the exciton (or tt particle) remains relatively large 
in the SC phase. The possibility of a critical point that 
is best described as corresponding to the regime exactly 
on the border between the bicritical and tetracritical be- 
havior was suggested in Ref. |107| (the projected SO (5) 
models discussed in that paper lead to such fine tuning 
for the effective theories). Other critical points, including 
a tricritinal-pne, have been suggested recently by Kivel- 



son et a/.E22l. We have argued here, instead, that many 
features of the experiments require the energy of the ex- 
citon to vanish at a quantum critical point describing a 
second order transition between the SC and SC+SDW 
phases; this appears when the finite temperature multi- 
critical point is tetra- critical (i.e. the four phases, SC, 
SDW, SC+SDW, and "normal" all meet at one finite 
temperature point) and has strongly broken equal-time 
SO (5) symmetry. We also note that AharonyE2j has re- 
cently shown, by an exact renormalization group analy- 
sis of fluctuations, that the finite temperature multicrit- 
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ical point has a 'decoupled' structure, which does indeed 
exhibit tetra-critical behavior. A finite coexistence re- 
gion between the superconducting and antiferromagnetic 
phases in the cup|rates has been also recently discussed 
by Martin et.ak 



iiZe have also mentioned the recent study of Kivelson et 
al.c^3 of a variety of finite temperature multicritical phase 



diagrams in three dimensions involving the SC and SDW 
order parameters. They pay particular attention to the 
possibility of a two-phase co-existence of SC and SDW 
order parameters, which should be distinguished from the 
homogenous SC+SDW phase we have discussed in this 
paper. In the presence of a finite field in the two-phase co- 
existence case, we would expect that the SC component 
has a H\n(l/H) term in its free energy, while the SDW 
component only has an analytic H 2 correction. Con- 
sequently, with increasing H, the fraction of the SDW 
component will grow at the expense of the SC compo- 
nent with an H\n(\/H) dependence. 

We mention that several other proposals for the ex- 
perimental consequences of the competition between the 
SC a nd SDW orders in the c uprates may be found in 
Refs. pTl|JTT2||n^jnip^||ll6| 



VIII. CONCLUSIONS 

The primary purpose of this paper has been a descrip- 
tion of the phase diagram in Fig y and of the static and 
dynamic properties of its low field phases. The point 
of departure of our work was the existence of a second- 
order quantum transition between the SC and SC+SDW 
phases in zero applied magnetic field (our methods can 
also be extended to weakly first-order transitions, but 
we did not discuss this here): we reviewed in Section || 
the early theoretical proposals and the experimental ev- 
idence in support of such a transition. In a non-zero 
field we found that this transition extended into a line 
of second-order transitions indicated by AM in Fig [|. 
This transition line approaches the H = axis with a 
vanishing derivative, which implied that relatively small 
fields could have a significant effect on the low energy 
spin fluctuation spectrum: this is our qualitative expla- 
nation for the field-induced enhancement pf low energy 
SDW correlations observed by Lake et alS3. Our anal- 
ysis also showed that the critical properties of the tran- 
sition in finite field were in all cases described by the 
familiar 0(3) symmetric ip 4 field theory: these have,.al- 
ready been described in some detail in earlier workuLi 
This mapping to the simple 0(3) continuum field theory 
occurs when the spin correlation length becomes larger 
than the vortex lattice spacing (as is always the case close 
enough to AM), and accounts for the fact that <Si a t pins 
the charge order fluctuations and so reduces the order pa- 
rameter to a real, 3-component vector. In principle, the 
Zeeman coupling to the 0(3) field theory modes should 
also be includeeLin the asymptotic critical region, but 
existing worktOtLZI has shown how to do this. We believe 



that experimental discovery of the critical field along the 
phase boundary AM is an exciting possibility for future 
investigations. Such a study should begin with a sam- 
ple with its s value slightly larger than s c ; application 
of a field should then allow tuning of the system across 
the quantum critical behavior associated with the AM 
phase boundary. The precise experimental control avail- 
able over the value of H should allow unprecedented ac- 
cess to an interesting, interacting quantum critical point 
in two dimensions. In the following subsection we dis- 
cuss a number of very recent experimental studies, and 
compare them to our results to the extent possible: we 
also mention proposals for future experiments. 



A. Implications for experiments 

So far, the most direct connection of our results with 
experiments is provided by neutron scattering measure- 
ments of the field dependence of the ordered moment in 
the SC+S|D.W| phase. Two such experiments have been 
performedljo in different but related compounds, and 
both show a reasonable fit to the predictedu H\nLL/H) 
dependence. The experiment of Khaykovich et alS3 ap- 
pears to be in a parameter regime similar to that of 
Fig [l^: there is an appreciable ordered moment at zero 
field, and the elastic scattering intensity roughly doubles 
in a field about a quarter of H c i- This is an important 
consistency check on our entire approach, as all numer- 
ical parameters in our computation had physically rea- 
sonable values. As is clear from Fig |l8|, the intensity of 
the satellite peaks associated with the reciprocal lattice 
vectors of the vortex lattice is quite small for these pa- 
rameters: this explains why such a satellite peak was not 
seen in the experiments even though they had the requi- 
site_wavevector resolution. The experiments of Lake et 
aZ.Efl are in a regime similar to that of Fig |l7]: they had 
quite a small moment at zero field, but this grew rapidly 
with field with a clear H\n(l/H) dependence. Again, 
as Fig 17 shows, the satellite vortex lattice peaks have 
a very small intensity, and this is presumably why they 
were not observed. This experimental sample appears to 
be rather close to s = s c , and we hope that a future ex- 
periment will move just past s c and study the transition 
across the AM phase boundary in Fig 

Our theoretical computations also suggest an approach 
by which the vortex (reciprocal) lattice may be detected 
in the spin fluctuation spectcijffl. While its influence on 



the elastic Bragg peaks 



was found to be very 



small in Figs [171 and 18, the spectra in Figs ffjlQ and E 



show a more significant influence in the inelastic neutron 
scattering cross-section. These plots may be viewed as 
the 'band structure' of the exciton moving in the vortex 
lattice, and the exciton dispersion shows clear features at 
the Bragg reflection planes in the reciprocal lattice of the 
vortex lattice. So we predict that a careful study of the 
inelastic neutron scattering spectrum may more easily 
yield evidence for the presence of the vortex lattice. 
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Next, we turn, to the recent STM measurements of 
Hoffman et alL3. These authors have observed sig- 
nals of charge order in the vortex lattice of BSCCO in 
the electron density of states at sub-gap energies. The 
charge order is at wavevectors K. cx — (7r/(2a),0) and 
K C y = (0, n/(2a)) (period of four lattice spacings), is 
peaked at the vortex cores, and extends about to a dis- 
tance which is about a quarter of the inter- vortex spacing. 
These measurements are most likely in the SC phase, 
where the SDW order is dynamically fluctuating. The 
nucleation of charge order by vortices in such a phase 
(but with the spins remaining dynamic) was predicted in 
Refs. |l0],[3l| Lattice scale theories^ of charge order in 
superconductors with preserved spin rotation invariance 
also found a substantial doping range of bond-centered 
charge order with a period of four lattice spacings, as 
did density matrix renormalization group studiea£3. The 
spatial extent of the envelope of this charge order in the 
SC phase has been computed in the present paper: the 
length scale in the observations isquitc similar to that in 
our computations in Figs [l5| and These computations 
were carried out for the same set of parameters (only the 
value of s — s c was changed to tune the doping level) 
used to obtain general quantitative consistency with the 
neutron scattering experiments above. The data of Hoff- 
man et al. seems rather similar to the result for f2(r) at 
point k in Fig [H| and the location of this point in the 
phase diagram of Fig ^| is very reasonable, given the op- 
timal doping of their sample and of their H value. This 
agreement suggests to us that the system studied by Hoff- 
man et al. has dynamic spin excitons, above a spin gap, 
which extend throughout the vortex lattice, as in Figs |l5| 
and |l6|; the charge order is then a signal of the pinning of 
these excitons by terms like those in 5i a t. An alternative 
model, in which the spin order was confined only to the 
region where charge order has been observed in STM, 
would have difficulty explaining the neutron scattering 
experiments: spin order so confined should yield easily 
observable satellite elastic Bragg peaks at the wavevec- 
tors of the reciprocal of the vortex lattice. 

Our computations also offer explanations for other fea- 
tures of the STM data which would be difficult to under- 
stand in terms of charge order nucleated independently 
in each vortex core: there is a noticeable correlation be- 
tween the phase and orientation of the charge order be- 
tween different vortices, which extends across the entire 
experimental sample. We believe this correlation is in- 
duced by the extended spin exciton states above the spin 
gap. Our model for the STM experiments can therefore 
be summarized as follows: the superflow in the vortex lat- 
tice reduces the energy of extended spin exciton states, 
and the sliding degree of freedom associated with spin 
density is then pinned by the vortex cores; this results 
in static CDW around each vortex, but the SDW order 
remains dynamic and gapped. A particular strength of 
our model is that it consistently explains the STM and 
neutron scattering experiments using the same set of pa- 
rameters. 



For the future, our theory suggests that neutron scat- 
tering and STM studies of SDW/CDW order should be 
carried out in systems where a uniform superflow has 
been induced directly by a current source, with no mag- 
netic field penetrating the sample. This will eliminate 
the vortex cores, but the superflow should still enhance 
the tendency for SDW/CDW order. Charge order can 
be pinned near impurities/defects of various kinds {e.g. 
dislocations, grain boundaries, surfaces), and so become 
visible to STM. 

We briefly comment on the high field phases (SDW 
and "Normal" ) in Fig ^, in which superconductivity is 
destroyed by the magnetic field. This regime may.be of 



relevance to the experiments of Boebinger et aLL_3. Dy- 
namic fluctuations of the superconducting order surely 
become important as we approach these phases, and so 
the theory of the present paper is not complete. Never- 
theless, given the nucleation of charge order near the vor- 
tex cores in ike SC phase (and its observation in the STM 
experimental^), it is natural to presume that this charge 
order survives into the "Normal" phase. The transport 
properties of the non-superconducting phases remain a 
very interesting topic for future research, but our naive 
expectation is that they are insulators. 

Another interesting type of experiments on supercon- 
ductors in thepiiprtex state has been, performed recently 
by Curro et. al&B and Mitrovic et. aZESE-D. They measured 
the local field dependence of the 17 spin-lattice relax- 
ation rate (1/Ti) and spin-echo decay rate (I/T2), this 
allowed them to deduce the rates as a function of position 
in the vortex lattice. Below we suggest how these experi- 
ments can be interpreted in our picture of the mixed state 
of the cuprates. The spin-lattice relaxation rate 1/Ti 
measures the rate at which nuclear spins are overturned 
as a result of interaction with electron spins. In. the BCS 
picture of vortices in a d-wave superconductorEj'Ef' ^ '' 
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this quantity is proportional to 7V(0) 2 and therefore in- 
creases dramatically close to the vortex cores due to sup- 
pression in the superconducting gap. On the other hand, 
as discussed in detail earlier in this paper, for the not too 
overdoped cuprates, charge density waves are nucleated 
around the vortex cores, which should lead to a suppres- 
sion in the local quasiparticle density of states, and hence 
1/Ti. This effect appears to have been observed in the 
experiments of Ref. |l6|. Another mechanism for the nu- 
clear spin relaxation is via the collective excitations of 
the electron system. In particular, the excitonic SDW 
excitations provide a large number of low energy S = 1 
excitations for flipping the nuclear spins. We suggest that 
a strong increase in the the high field part of 1 jT\ (cor- 
responding to the vortex cores) with increasing magnetic 
field in the experiments of Mitrovic et.al&B reflects the 
growth of the SDW correlations and the corresponding 
increase in the excitonic susceptibility. It would be in- 
teresting to study this enhancement quantitatively and 
compare it with the H In (1/77) behavior observed in 
neutron scattering experiments and derived theoretically 
in this work. We mention that the non-two-sublattice 
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SDW makes this mechanism more effective for relaxing 
the 17 nuclear spins, in contrast to the (ir, tt) electron 
magnetism which leads to a magnetic field on the oxy- 
gen sites only through the Dzyaloshinskii-Moriya inter- 
action and weak ferromagnetism. The echo decay rate 
I/T2 is related to the inhomogeneity of the local mag- 
netic fields. The appearance of the local SDW order (or 
sufficiently slow fluctuations) should therefore contribute 
to the increase in I/T2. The SDW enhancement is rel- 
atively stronger around the vortex cores, which should 
give rise to the enhancement in 1/T 2 in this region; this 
agrees with the experimental observations in Ref. [l5|. The 
analysis of our paper suggests that the difference in 1/T 2 
will not become very large upon approaching the SC to 
SC+SDW boundary, as the SDW excitations become ex- 
tended close to this phase boundary. As the magnetic 
field is increased, the SDW fluctuations should become 
more pronounced, so we expect that I/T2 will increase for 
all values of the local field. By contrast, in the BCS the- 
ory, one would expect that I/T2 decreases with increasing 
magnetic field, since the field becomes more uniform. 
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APPENDIX A: DZYALOSHINSKII-MORIYA 
INTERACTION 



An orthorombic distortion of La2-sSr$CuOi results 
in the Dzyaloshinskii-Moriya (DM) interaction for the 
Cu spins, 



n 



DM 



= \J2(~yd-s t x s t+5 , 



(Al) 



i..S 



where the sum over S extends over all the nearest neigh- 
bors of site i, and d is a unit vector in the direction^ 
the orthorombic a axis (i.e. a diagonal (1, 1) direction)E3. 
In this appendix, we study the effect of the DM interac- 
tion on the non-two sublattice SDW ,and for simplic- 
ity we consider a SDW at one wavevector only. The 
Hamiltonian (Al) mixes wavevectors q and Q + q, wh ere 
Q = (n/a, 7r /a). In this case we need to modify (1.2) to 



S(t,t) = Re 



r C(r, r) + e^+^M^r, r)l(A2) 



Straightforward algebra shows that the contribution of 
the DM interaction to the action is 

S DM + S M = J d 2 dr { X(K S ) d.[i,x M* + c.c] 



2x 



} 



(A3) 



where A(K S ) = 2A(cos(K s:E a a ;) + cos(K S j,aj,)) and the 
last term comes from the fact that spin fluctuations M x 
are massive. We can now integrate M out, and find the 
anisotropy term for the SDW order parameter 



A 2 (K S ) X 



d 2 dr\ $ x xd\ 2 



(A4) 



From (A4) we see that the DM interaction favors the 
collinear SDW, with direction of <5 perpendicular to d, 
i.e. along the orthorombic b axis (direction of the SDW 
ordering is always in the CuO plane). We also expect 
that the anisotropy becomes weaker with increasing dop- 
ing due to a decrease of A(K,jL However, the typical 
scale for the anisotropy is smallEZI, and so we expect that 
the quartic it2|<& 2 | 2 term in (1.8) plays a dominant in se- 



lecting the collinear SDW at low temperatures. We note 
that Sdm is quadratic, so it will favor the collinear SDW 
fluctuations even above the transition temperature. 



APPENDIX B: ZEEMAN COUPLING TO THE 
MAGNETIC FIELD 

This appendix briefly discusses the effect of the Zee- 
man coupling to the magnetic field on the action iS$ in 
(1.8) for the SDW fluctuations. We will see that the ef- 
fects are weaker than those considered in the body of the 
paper, especially near the critical point A at s — s c in 
zero field (see Fig ||). 

As reviewed in Ref. |121| , in systems without an over- 
damped particle-hole continuum of spin excitations (as 
is the case here hear the ordering momenta K SX)2/ ), we 
can deduce the coupling to the external field using simple 
gauge invariance arguments. In particular, the primary 
consequence of the external field is to rotate the spins 
uniformly about the field axis, and this can be accounted 
for by the following replacement to all temporal gradient 
terms: 



(Bl) 



and similarly for <f> yQ . Here H a is the three vector in spin 
space representing the external field. The resulting <S$ is 
closelM *&Lated to models that have been studied in some 
detailQ'Lill in the context of double layer quantum hall 
systems. From this work, we can deduce the phase dia- 
gram sketched in Fig |3^. The most important property 
of this phase diagram is that zero field phase transition 
at s = s c moves to finite field as H ~ (g — s c ) zv where 
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spin-space plane 
orthogonal to the field 
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FIG. 32: Phase diagram of 5$ in (1.5) including the Zeeman 
coupling in (|Bl[). The point A is the same as the correspond- 
ing point in Fig ^| The central argument of Appendix |§| is 
that it requires a much larger field for s > s c near A to induce 
SDW order above, than in Fig [| 



the exponent zv = 1/2 in mean-field theory. Fluctua- 
tion corrections will slightly increase this value, but the 
present critical field will nevertheless remain larger than 
the field in (3.9) associated with the corrections arising 
from the superflow. In particular, the phase boundary in 
Fig [32] approaches the H = line with an infinite slope. 
Consequently, the Zeeman shift is subdominant to the 
stronger effects discussed in the body of the paper. 



APPENDIX C: MICROSCOPIC THEORY FOR 
COUPLING BETWEEN SC AND SDW ORDER 
PARAMETERS 

In this Appendix we discuss the microscopic origin of 
the effective interaction k between the SC and SDW order 



parameters in (1.10). We will argue that repulsive n > 
is a remarkable property of doped Mott insulators, but 
weakly interacting electron systems quite possibly have 
k < 0. 

We start by considering a weakly coupled Fermi liquid 
of electrons Ci a moving on the sites, i, of a square lattice 
which is close to superconducting and commensurate an- 
tiferromagnetic instabilities: 

Z = J V^Vce s[c] 

s[c] = dT^24d T a-n[c]^ 

H[c] = efccj^Cfcg + Himp 

ha 

Al 



k 
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k/3 
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-co, -p 
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FIG. 33: Diagrams that contribute to the effective interac- 
tion between the superconducting and antiferromagnetic or- 
ders for Fermi liquids. Solid lines correspond to the quasi- 
particle propagators, zigzag lines correspond to the d-wave 
superconducting order parameter and contribute a factor d v ; 
wavy line describe the SDW, and a dashed line describes the 
static disorder potential. 



conducting d-wave order parameter, and we assume a 
nearest-neighbor tight binding dispersion of the electrons 
ej; = — 2i(cos /c^ +cos fcj,) — /1, Q = (7r,7r), and everywhere 
in this section momentum integrals go over the first Bril- 
louin zone. 7Yi mp describes the static potential of the 
impurities which gives rise to a finite quasiparticle life- 
time 



1 



27rn imp A^(0)V- 2 , 



(C2) 



where N(0) is the density of states on the fermi level and 
V is the impurity potential. 

Assuming that A and $ are small we can integrate 
out the fermions and obtain 



-0F 



F = F GL [Ao} + F AF [$} + K \Ao\ 2 $ 2 



(C3) 



The diagrammatic representations of the terms that con- 
tribute to k are shown on Fig. 33. Solid lines correspond 
to the quasiparticle propagators 



G(p,u> n ) = (iuj n - e p - i/(2r)sgn(cj)), 



(C4) 



Here = Ao(cosfc x — cosfcj,)/2 = Aodk is the super- 



zigzag lines correspond to the c?-wave superconducting 
order parameter and contribute a factor d p ; wavy line 
describe the SDW, and dashed line describes the static 



31 



disorder potential. We have 
1 f d 2 p 



(a) 
(b) 

(c) = 

with 

L(aj v 



pUp+Q 2^ 



PJ (2^) 2 
x G(u n ,p + Q)G(-cj n ,-p-Q) 

l r d 2 P 

PJ (2tt) 2 



G(cj n ,p)G(-w n ,-p) 



G(w„,p)G(-u n ,-p) 



x G(w„,p + Q)G(-Lj n , -p - Q) 
V 2 ^y^L{cu n )M(Lo n ) 



(C5) 



d 2 p 

d 2 G(uj n ,p)G{-u n , -p)G(uj n ,p) 

d 2 p 

G(u n ,p)G(-uj n , ~p)G(uj n ,p + Q) 

(2ir) 2 

(C6) 



It is useful to note that if we define the static spin 
susceptibility at momentum Q in the superconducting 
state 



x(Q) 



d 2 p 

+ F(u n>P )F(ui n ,p + Q)}, 



{G sc {Un,p)G S c(Un,P + Q) 

(C7) 



be directly applied. Instead, we have to turn to alterna- 
tive strong coupling approaches, in which the existence 
of the Mott insulator is built in at the outset. Such ap- 
proaches have been discussed recently, and these are ex- 
pressed in terms of collective degrees of freedom which are 
natural in the vicinity of of Mott insulator. Electron spin 
singlet states, spin one magnons, Cooper pairs of holes, 
and fermionic quasiparticles are introduced as individual 
excitations, and interactions betweenJiem are obtained 
from the microscopic t-J HamiltonianEJlHj (phenomeno- 
logical models of just Xhe bosonic degrees of freedom have 
also been considered!!^!!) . All these papers find strong 
repulsion between magnon and hole pair states, arising 
from the constraint on the allowed Hilbert space. The 
origin of this repulsion therefore lies in the short distance, 
lattice-scale physics of allowed low energy states near a 
Mott insulator, rather than effects near the Fermi sur- 
face in the weak-coupling analysis discussed earlier. The 
Cooper pair-magnon repulsion immediately implies re- 
pulsion between the superconducting and antiferromag- 
netic orders, since the superconducting and antiferromag- 
netic phases correspond to the condensates of the corre- 
sponding particles. As an example, see Fig 2 in Ref. |5|: 
the pairing amplitude is weak in the region with mag- 
netic order, but rises rapidly once the magnetic order is 
suppressed. 



with the Green's functions in_the superconducting state 
defined in the usual mannerliH 2 ], then 



d 2 x(Q) 
dA*dA 



A o =0 



(C8) 



which agrees with ( |C5| ) and (p6|). 

In the limit /it 3> 1 the main contribution to k comes 
from the diagram (a) and we find for T — > 



K = -7T<4>iV(0)~£ 



1 



^ (w, 

r n>0 v ? 



+ l/2r)(/i 2 + K + l/2r) 2 ) 
(C9) 



It is important to note that in deriving the expression 
( |C9| ) we relied on the fact that we have a c£-wave super- 
conductor with d p+ Q — —dp and took the average value 
of (d 2 ) on the Fermi surface to be 1. Hence, such Fermi 
liquids on the square lattice have an effective "attraction" 
between the antiferromagnetic and superconducting or- 
ders, which can be traced back to the enhancement of 



the antiferromagnetic susceptibility (C7) in the e?-wave 
superconducting state. 

We have so far examined the interplay between SC and 
SDW orders near the boundary of their instability to a 
weakly interacting Fermi liquid. Now let us turn to the 
same interplay, but in the vicinity of a Mott insulator. 
Strong interactions are required to produce the Mott in- 
sulator, and so the perturbative approach of ( pl| ) cannot 



APPENDIX D: RENORMALIZATION GROUP 
ANALYSIS OF COMPLEX VECTOR FIELDS 



This appendix will briefly review existing theoretical 
results for the critical properties of field theories which 
are similar to 5$, but simpler. The analysis of the full 
5$ theory will be addressed in future work. 

The simplification made here is to consider a field 
theory with only one complex vector field, $ a , with 
a = 1 . . . m; the original model has two such fields $ XQ 
and $y : a- For the case of only one such field, we can 
always rescale x and y co-ordinates to make all velocities 
unity; then in d space dimensions we are interested in the 
field theory with action 



Sr. = 



d d rdr 



|d T $ a | 2 + |V r $ Q | 2 + S |$ a | 



y|4> Q | 4 + yl<^l 



(Dl) 



This theory has upper critical dimension d = 3, and 
can be studied in an expansion in e = 3 — d. Renor- 
malization group equations for the quartic terms were 
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obtained to 0(e 2 ) by D. Jones et. all 



BB0 



111 



2Ki 



- (3m + 7)u\ + 22u\u 2 



du2 



+ (5m + 24)mu\ + 4(m + 2)u\ 



= eu2 — Kd[mu\ + &U1U2] 

- 2K I (m - A)u\ - 2(5 + 3m> 2 ui 
- i(5m + 41)u 2 u 2 (D2) 



where K d = 2- d+1 /[ir d / 2 T{d/2)\. These flow equations 
always have two unstable fixed points: the Gaussian 
2 = and the isotropic 0(2m) Heisenberg 



point u* — 
fixed point 



e 

K d m 
= 



Q(e 2 



(D3) 



For sufficiently large or small m there may also be two 
other fixed points 



u 



QK d 

v* = —B 
K d 



B m [3m 2 - 12m + 144 qp 3mi?^ /2 ] e + 0(e 2 



m 



12±3R^ 2 }e + 0(e 2 ) (D4) 



where i?" 1 = m 3 + 4m 2 — 24m + 144 and R m = m 2 — 
24m + 48. The last two fixed points are absent in the 
case of m = 3. We note, however, that for m = 2 and 
in the large m limit a stable fixed point (the so-called 
chiral fixed point, see Ref. |37j) is possible for U2 > 0, so 
it may control the transition to the spiral order. When 
U2 < 0, the system always flows towards strong coupling 
U2 — * —00, so we expect the transition to collinear order 
to be weakly first-order. 



APPENDIX E: NUMERICAL SOLUTION IN THE 
SC PHASE 



We wi 



BrandtEi 



use the methods and notation described in 



First, assume we know Vfr(r). Write its Fourier ex- 
pansion in the form 



G r 



(El) 



G 



where da = gLg are both real, and G are the recipro- 
cal lattice vectors of the triangular vortex lattice. Un- 
like the convention followed by Brandt, the sum over 
G always includes G = 0, unless stated otherwise ex- 



plicitly. In Brandt's notation, (El) can be inverted by 



d G = (V,f/(r) cos(G • r)), where the angular bracket de- 
notes a spatial average. Because of the symmetry we can 
work on only half a unit cell of the vortex lattice, and 
for simplicity we choose the half unit cell to be the one 
plotted in Fig. [l^. 

To obtain Gh, we want all the eigenvalues and eigen- 
functions of the Schrodinger equation ( 4. 20] ). As in the 
usual Bloch theory, these are labeled by a wavevector k in 
the first Brillouin zone, and a band index fi. The explicit 
form of these are 



s ^k(r) 



e 4k r ^ 



CuG 



G 



(k) 



AGr 



(E2) 



where Au is the area of the unit cell, and the c M G(k) are 
normalized so that 



(E3) 



If we choose M values of G (also as in Brandt), then fi = 
1 . . . M , and the c^g (k) are the orthonormal eigenvalues 
of the M x M matrix AlG,G'(k) where 

^G,G'(k)c^(k) = £ 2 (k)c MG (k) 

G' 

A4g,g< (k) = (k + G) 2 ,5g,g' + d G - G , (E4) 

After this diagonalization we obtain the Fourier compo- 
nents of (El) as 



do = s — s 

:[ 



G 



^pth(ir G (k)/(2rn 

k,G 



Nu 
N k A u 

1 



A 2 
^0 



,(E5) 



and for G ^ 
d G = 



Nu 



kq g 



CuG' 



2 NkAu J-i '' 



%(G'+G)(k) 



-C M (G'-G)(k) 



3th(£ M (k)/(2r)) 



(E6) 



where the sum over k is over Nk points which average 
over the first Brillouin zone. Also note that c^ G (k) = 
c M 7^(G)(7?-(k)) where 1Z denotes a rotation by 7r/3. This 
can be used to cut the numb er o f k points in 1/6. 

The iteration of ( fE5| ) and ( |E6|) will produce the solu- 



tion to (4.17) for a given 



The next step is to solve ( |4. 18) ) , given V#(r) in (|El|). 
This is done just as in Brandt. His Eqn (9) is replaced 
by 



4uT 



(-V 2 + 2)lu = 2 



(Vh~A 2 ) 



—) 

4uT J 



IV'ol 2 - 



c^Q 2 - g }(E7) 
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and a corresponding change to Brandt's Eqn (11) 
new form of Brandt's Eqn (12) is 



The and 



4uT 



a G x 
l^ol 2 



4uT 



(Vir-Ag) 



(w 2 )(1-k 2 /(4uT)) 



(E8) 



After determining lo from abo ve, w e use this result 
to obtain new Vg( r) by solvi ng (4.17 ), and so on. By 
iteration of Eqn. (4.17) and (4.18|), we will be able to 



have the final solution to both of them. 

Note that in order to get our numerical results, we did 
use a finite momentum cutoff. However, the equations 
have been designed to be cutoff independent and we did 
find that the Fourier components of ^ij(r) and Vff(r) 
decreases rapidly upon going to higher momenta. 



APPENDIX F: SPIN ORDERING PHASE 
BOUNDARY NEAR M 



Her e we discuss the analytical solution of ( 4.17 ) and 
(4.18) in the vicinity of the multi-critical point M in Fig||, 
with the aim of determining the location of the AM phase 
boundary in its vicinity. Analytical progress is possi- 
ble because the amplitude of the superconducting order, 
\iPh{y)\ 2 i is small in this region. Our analysis will show 
that in this region AM behaves as H = 1 — q(k — s + s c ), 
where g is a numerical const ant. The earlier full numer- 
ical solution in Section IV C led to the estimate pa 1.2, 



and we shall find a consistent result here. 



In addition to (El), we use the Fourier expansions 



ducPk 



P- + k 2 + A 2 , 

G 



(Fl) 



b G = - 



IG 

4tt 2 



d 2 k 



(k + G) 



coth(Vfc 2 + d /2T) 



coth( ^/(k + G) T +do/2T) 



G) 



IG 



2^{k- 



d G . „-,, ,2 



- d Q 
at T = 



2Vfc 2 + d 



G ^ (F5) 



Now we can solve (F3,F4,Fj|) for the d G in terms of the 
aG- 



Finally, we need to determine the a G by solving ( 4.18 ). 
This can be done with the realization that for small ipn, 
the functional form of the superconducting order parame- 
ter can be assumed to be equal to the Abrikosov solution. 
So we assume 



aoa G 



G ^ 



(F6) 



where Bq is given in (8) of Brandt. Now, it remains to 
obtain a single additional equa tion t o determine cio- This 
we determine by multiplying ( 4.18 ) by ^lr( r ) an d aver- 
aging over all space. Using the property of the Abrikosov 
solution for ?/>#(r), we obtain 




«g«-g - \i>o\ 2 a 



Ha = 



(F7) 



Eqns (F3|,F4,F5 F6,F7) are now simple equations that 
can be easily solved to obtain all the Fourier coefficients. 
The line AM corresponds to do = 0. Our analytical result 
of the slope of AM near M point is g w 1.1, which is in 
acceptable agreement with that obtained from the full 
numerical solution. 



|<Mr)| 2 =^> G( 



tGr 



(F2) 



G 



Note that this notation for a G is slightly different from 
that above and in Brandt. 



Then (4.17) becomes 

d Q = A 2 + K(a - \ifj Q \ 2 ) + 2Nub Q 
d G = na G + 2Nub G ; G ^ (F3) 



Second, we can solve ( |4.q ) by a Feynman graph expan- 
sion in c?g^o- This yields 



, kdk 
= 1 ^ 
1 



coth(VF + 5o/2T) 



2^/k 2 T 
<%- A 

47T 



A 2 , 



2VF + 5o 



-0{d 2 G ^) atT = (F4) 



APPENDIX G: NUMERICAL SOLUTION IN THE 
SC+SDW PHASE 

H er e we will des cribe the solution of the equations 
(0), Q and (PI for the unknowns V H (r), ip H (r) and 
Tijj(r). First, as (5.4) is linear in n#(r), it is convenient 
to rescale 



n H (r) -> n H (r)/V2Nu, 
and these equations become: 

V h (t) =s-s c + n{\^ H {v)\ 2 - 1) + n 2 H (r) 

r x -> „ , f dujd 2 k 1 

+ 2Nu[T^G H (r,r,. n) - ^ ^ 

(-V 2 +V ff (r))n H (r) = 0, 



(Gl) 



,(G2) 



(G3) 
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4uT 



l-A,)(hMr)| 2 -l)+ '' 



;(Vh(t) - S + S e ) 



Nu 



4uT 

(V r -zA) 2 ]^(r) =0(G4) 



We use t wo-step iter ation to self-consistently solve the 
equations (G2, G3|, |G4| ). The first step consists of solving 
Eqn. (G2) and flG3| )"7"and the second step is solving (|G4|). 

For the first step, we use a four-substep iteration. 
First, define and calculate 



+2Nu\^T^Gh{ 



s-s c + k{\iPh[v)\-1) 
cLucPk 



r, r,w„ . - 



1 



8tt 3 uj 2 + k 2 



Second, define and calculate the inverse of operator 

A = -V 2 +r] H (r). 
Third, calculate njy(r) which satisfies 

n H (r) = -A^n^v). 
Last, calculate Vir(r) using 

V H (r)=VH(r)+n 2 H (r). 



.(G5) 

(G6) 
(G7) 
(G8) 



Choo sing p rope r i nitia l valu e for Vjj(r) and n#(r) and 
iterate ( U5h , ( |G6| ), (G7), (|G8|) will produce the solution 
to both (|G|) and (|G3^ 

In practice, the above steps are performed in momen- 
tum space. If we let 









G 




= E /<^ Gr 




G 


1//W 


= E^ iGr 



(G9) 



G 



where G are the reciprocal lattice vectors of the vortex 
lattice (note that fc dif fers slightly fro m /g in ( 5.10| ) 
because of the rescaling (|GlD), then (G5 - |G8|) become 



flo = s 
Nu 



S c + k(^«G - 1) 



G 



N k A 



u k+G#0 



coth(g G (k)/(2T)) 



Ik + GI 



(G10) 



9g = -kog 



2N k A u 

C M (G'-G)(k) 



E E 



k,G' M ,B M (k)#0 

coth(£ M (k)/(2T)) 



C i u(G'+G)(k) 

(Gil) 



-^GG' — G 2 Sg,G' + .9G-G', 



(G12) 



/ G - -E^GG'KWcos(G'.r)), (G13) 

G' 

d G = g G + (nl(r) C os(G' -r)). (G14) 
Note that in the substcp (G13) the equation is solved 
by another smaller iteration. 

The second step is very similar to the case with no 
magnetic order as in Appendix [e| The equation (G4) 
can be solved by a two-substep iteration of the following 
equations 



(-V 2 + 2) u; = 2 (1 + G(r))w 



(G15) 



where 



G(r) 



2 \ 



1 4uT J 4uT 



(V ff (r)-s + s c ), (G16) 



and 



(ujG - ^g 2 - ff) 
(w 2 )(1-k 2 /(4uT))' 



(G17) 



From the iteration results we are able to determine 
Vff(r), rifl(r) and V'ff(r). 
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